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1-0) Preface 

The development of the Space Shuttle will provide a unique vehicle 
from which the global distribution of aerosols may be measured. The Shuttle 
missions will occur at a time when the interests of the meteorological 
community are concentrating on a better understanding of the earth's climate. 
Desert aerosols are the most abundant atmospheric aerosols and it appears 
that they will significantly affect the climate. A Shuttle-borne lidar 
system can provide basic data about aerosol distributions for forthcoming 
models of climate. 

Chapter 1 of this report summarizes our present knowledge of the physical 
characteristics of desert aerosols and the absorption characteristics of atmo- 
spheric gas. Chapter 2 presents an analysis of radiative heating computations. 
Computational simplicity and adequate accuracy are fundamental criteria which 
any computational scheme must satisfy in order to be useful as part of a com- 
prehensive climatic model. Sample computations of the radiative heating effects 
produced by desert aerosols are also presented in Chapter 2. General circulation 
models are insensitive to changes in heating rates less than 0.5°K/day; this 
ultimately determines che accuracy with which the aerosol distribution must be 
determined by the Shuttle-borne lidar. These considerations are presented 
in Chapter 3- 

A high spectral resolution lidar can provide profiles of the aerosol 
extinction coefficient. Tables 3-1 and 3-2 establish criteria for height 
resolution and the accuracy with which the lidar must provide aerosol extinction 
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i- coefficients. The analysis of Chapter 4 presents the characteristics of a I 

f: i 

?? Shuttle-borne lidar that will satisfy the criteria established in Chapter 3. , 

X j 

h Chapter 4 also presents some laboratory studies which seek to identify j 

t: 

k schemes that might permit the implementation of a high spectral resolution s 

h $ 

& s 

r lidar system. 
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1-1) Introduction 


Climatic variability is a major hazard which afflicts agricultural 
production. Much attention has been devoted to the influence of atmospheric 
constituents on the heat budget and its possible long-term impact on climatic 
trends, see e.g. SMIC report V 197 1 ) . The distribution of radiatively active 
constituents such as aerosols and clouds may produce changes in the distri- 
bution of surface temperatures (duration of growing season) or precipitation 
(drought or floods) before longer range trends such as ice ages become sig- 
nificant. General circulation patterns are influenced by the vertical sta- 
bility of the at':c sphere. The temperature profiles wMch determine this 
stability are affected by the vertical distributions of radiatively active 
constituents which may alter radiative heating or cooling rates. 

Global measurements of the vertical distributions of clouds and aerosols 
will be required as input for models needed to study the influence of aerosols 
and clouds on climatic variability. Lidar can be used to measure the vertical 
distribution of aerosols and clouds; and the development of the space shuttle 
provides a vehicle which may be used to transport a lidar system over much of 
the earth. 

Cloud top heights are easily determined by lidar. The optical thickness 
of cirrus clouds can also be readily determined by lidar. While these quanti- 
ties are of considerable meteorological interest in their own right, we will 
confine our attention to aerosols and their heating effect on the atmosphere. 
The reasons for this choice are that desert aerosols are more wide spread 
than anthroprogenic aerosols, they are found nearer to the equator where 
they may interact more effectively with solar radiation, they frequently 
appear in cloud free regions so that their radiative effects are not obscured 
by clouds, and they are more difficult to measure with a lidar than clouds . 


Hiw«» i~ — . : 
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1-2) The Desert Aerosol 

1-2-a) Geographical Distribution and Climatic Role of Deaert Aerosols 

The desert aerosol is probably the most vide-spread natural aerosol. 
Figure 1-1 shows the approximate global distribution of desert aerosol sources. 
The largest source is the long band of deserts and steppes stretching from 
West Africa east through the Sahara, the Sahel, the Lybian and Egyptian 
deserts, the Sudan, Saudi Arabia, Jordan, Israel, Syria, Iraq. Iran, 

i 

! ' Pakistan, and India. The second source is located in the extensive Siberian 

■ • 1 

and Chinese deserts. Other sources, are located in South Africa, South 

! 

i l i 

America, Australia and the Western U.S. All together, deserts comprise about 

i! 1 i 

eight percent of the Earth s surface or close to one third of the total land 

(' 

surface. The total output of "long-lived" dust from this source is at least 

i 

120 million metric tons annually according to Carlson and Prospero, (1972) 

\ 

and Joseph et al., (1973) and it may reach more than four times this value, 
see Goldberg, (1971). An upper bound to the global atmospheric aerosol load 
of the order of 2000 million metric tons/year was cited in the SMIC report, 
(1971). 


There is evidence that overcultivation of marginal lands promotes the 
development of deserts. As population pressure places more marginal land 


under cultivation, the aerosol density in the atmosphere may be expected to 
increase. 

One reason for the climatic importance of desert aerosol is that small 


changes in climate or in man’s activities may cause considerable dis- 
placement of arid zone boundaries. The more arid the climate thus becomes, 
the more desert aerosol will be introduced on a global scale into the atmo- 
sphere. There might thus be an initial positive feedback mechanism between 
desertification and a change in ra.-.'.-.tJon balance due to desert aerosol. 


- ^ 


. \ 




Tig. 1-1 Global distribution of desert aerosol sources and regions over which the atmosphere 
frequently contains aerosols. From The Times Atlas of the World, London, Times Nevspapers Ltd. 
(1968) pp 272. 
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Heavy desert dust deposition far away from its source has been reported 
extensively. North African dust has been reported in Northern Europe by 
Morikofer, (1941) and Stevenson, (1969), in the West Indies by Prospero, 
(1968); and Carlson and Prospero, (x971). Desert aerosols have been shown to 
exist in the atmosphere at low and sub-tropical latitudes by Peterson, (1968) 
Bryson and Wendland, (1970); Volz, (1970); Volz and Sheehan, (1971); SMIC 
(1971); Joseph and Kanes, (1971); Joseph et al., (1973); and Sargent and 
Beckman, (197. 1 ,) . 

The appearance of desert dust in large quantities over Japan and the 
Western and Central Pacific is correlated in spring with the northward migra- 
tion of the 0°C isotherm in the Central Asian deserts. T. Kitaoka, (1971), 
indicates that the turbidity over Japan may increase by 30CZ when this dust 
appears. The wind-borne desert aerosol from Asia has also been found in 
geologically and pedological ly important amounts on the islands of Hawaii 
by Jackson et al., (1971). Heavy dust falls reported by ships in the Pacific 
up to 3000 km from the source, have been cited by Ing (1972). 

The horizontal distribution of aerosol in large scale air masses is 
dense enough so that its presence is easily detectable in broadband satellite 
photographs obtained in reflected sunlight. 

1-2-b) Vertical Distribution of Desert Aerosols 

Desert aerosols are very dense over source regions. Surface concen- 
trations are very frequently of the ord.rt of 1000 Jjg/m^ according to 

Peterson, (1968). The average surface mass concentration in the desert 

3 

aerosol layer is between 30-300 pg/m away from the source, see Carlson 
and Prospero, (1972), Reynolds et al., (1973) and Joseph et al. , (1973). 

These mass concentrations are comparable to those found where man-made 
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pollution is severe. This is one reason why natural desert aerosols nay be 
more significant to the global energy balance than man-made aerosols. 

The presence of desert aerosols is clearly evident in visibility measure- 
ments. In very dense desert haze (not a sand storm) the visibility may be less 
than one ks. In the Middle East area, days with visibility less than about 
10 km have a frequency of occurrence of about 50-200 per year according to 

Manes and Guetta, (1973), and Lentz and Hoidale, (1974). 

. \ 

| Prospero, (1968), and Carlson and Prospero, (1972) found that desert 
aerosols are usually transported away from their source in an lsentropic layer. 

;i - ! 

Over the Atlantic Ocean this layer is usually found between 800 and 500 mb. 

! i , ; 

Below this layer, the air contains very little desert aerosol, Junge, (1972). 
The same stable layers in which the desert aerosol is advected were also 
found over the Middle East by Joseph, Ashbell and Eviatar, (1961). However, 
the aerosol has been found well-mixed to high altitudes over the source region 
by Peterson (1968), and Peterson and Bryson, (1968). 

1-2-c) Physical Properties of Desert Aerosols 

A study of the effect of a model aerosol, with properties similar to 
a desert aerosol, on total solar and infra-red radiative fluxes, flux- 
divergence and radiative-convective temperature profiles, has been published 
recently by Wang and Domoto, (1974). They found that the absorption of the 
aerosol was equal to that of the water vapor in the solar spectral region. 

Over regions of high albedo (e.g. deserts), an increase in the aerosol concen- 
tration will heat the earth-atmosphere system. Furthermore, convection is 
suppressed by the presence of aerosol and in the case of a "dense" aerosol 
(visibility ^ 5 km), a temperature inversion lsyor is formed. Net heating 
rates in the lowest layers are changed by about 2.0°K day 
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An experimental study by Reynolds et al., (1973) shows that the presence 
of a layer of desert aerosol aloft In the tropical oceanic atmosphere may 
shift the net (solar + thermal Infra red) radiation heating rate curve downward 
in altitude and Increase the heating at the altitude where the concentration la 
a maximum. 

The stabilizing effect of solar absorption by the atmospheric aerosol 
and its greenhouse effect in the thermal infra-red may play a significant role 
in maintaining temperature inversions over deserts and it may affect convection 
over oceans. The peculiar height distribution of the desert aerosol over the 
tropical oceans may increase the static stability enough to inhibit the 
rainfall during day 1* gin. hours. 

Joseph, (1970), (1971), Stowe, (1971), and Jacobow^tz and Coulson (1973) 
have shown that the presence of aerosol over a surface will also significantly 
affect the remote determination of the surface temperature from a satellite. 

It is thus extremely important to study the vertical distribution of 
the desert aerosol on a global basis. 

1-2-d) Composition of the Desert Aerosol 

The composition of desert aerosols has been investigated by Junge, (1972), 
Goldberg, (1972), Carlson and Prospero, (1972), Kondratiev, (1972), Volz, 

(1971), Peterson (1968), Lentz and Hoidale, (1974), Llndberg and Laude, (19/5\ 
and Eullrich et al. (1973). 

The largest fraction of the desert dust is a variety of silicate clays, 
together with admixtures of calcite, quartz, sodium nitrate and organic matter. 
Organ natter or iron compounds may be necessary to explain the relatively high 
measured absorptivity of the desert aerosol for visible solar radiation found 
by Joseph and Wolfson (1975), and Kondratiev (1972); (1974). 
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The composition by weight summarized in Table 1-1 is based primarily on 


work by Lentz and Hoidale, (1974). 


Table 1-1 Composition of Desert Aerosol 

Montmorillonite 

352 

Kaolinite 

20 % 

Illite 

20% 

Calcite 

10% 

Iron Compounds (Limonite), Organic Compouncs 

5% 

1 Quartz 

5% 

1 

■I Sodium Nitrate 

;l f, 

!•' ‘ S 

5% 


| The sllicaceous clays that comprise the bulk of desert aerosol and 

\ 

possibly organic material that adheres to the clay particles, found by 
Gillette et al., (L974), are excellent natural clcud seeding agents. This 
fact may lead to a significant effect on cloud and rain formation and cloud 
optical properties. This property of the clays may lead to t profound 
diurnal change of the aerosol over deserts. During daytime, when the relative 
humidity becomes quite low, aeroso.s may be dry and non-spherical. Such 
aerosols absorb both solar radiation and thermal radiation. At night the 
relative humidity often approaches 100 percent; the aerosols may rapidly 
change their optical properties, their shape and their size distribution as 
water is adsorbed onto the aerosols. Such aerosols will absorb less in the 
visible in the early morning hours till the water-vapor evaporates. Some 
experimental Indications of such a phenomenon have been reported by 
Vittori et al., (1974) for aerosols in Italy. 
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Lindberg, (1974) found that the hydroxyl bonds In clays have absorption 
bands in the near infra-red and these may simulate the radiative effect of 
water vapor. The absorption of solar radiation by desert aerosols is about 
equal in magnitude to that of water-vapor according to Reynolds et al. t 
(1973); and Kondratiev, (1972,1973), Kondratiev et al., (1972), (1974). 

1-2-e) Refractive Index of Desert Aerosol 

(1) The Solar Spectral Region: .4 <im < 1 < 1.1 pm. 

Most measurements and theoretical evaluations agree that the real part of 
the refractive Index of "dry" (ambient relative humidity ^ 60%) desert aerosol 
is of the order 1.46 < n < 1.54. This is the range of values appropriate 
for both amorphous and crystallite quartz according to Day et al., (197&), 
and for siliceous clays like kaolinite according to Steyer et al., (1974). 

It is too low for iron or organic compounds which comprise at most 5% of 
desert aerosols. We suggest the value of n = 1.54 in the spectral range 
.4 pm < A < 1.1 pm. 

There is much less agreement on the value of the imaginary part of the 
refractive index, k, see Fig. 1-2. The highest values are those deemed 
appropriate for limonite, see Kondratiev et al., (1972),(1974) . The lowest 
values are those of Lindberg and Laude, (1975), based on reflectivity measure- 
ments of a desert aerosol diluted in a powder matrix. Brinkworth (1972) 
suggests that such results may be too low by a factor of about 2. The values 
used by Mc.Clatchey et al. (1970) are based on Volz's (1957) data. Joseph and 
Wolfson (1975) deduced a mean value for the spectral region .3 Pm < A < 2.5 pm 
of k .03 ± .02 on the basis of solar flux measurements. The curve 3K In 
Figure 1-2 shows the Lindberg and Laude curve shifted upward by a factor of 
3 to coincide with that reported by Kondratiev. The range of uncertainties 
for k thus seems to lie between that of limonite and the "desert dust" curve 
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of Lindberg and Laude (1974), in which there is no limonite. It is Impossible. 
to fix the value any better. The curves presented by Kondratiev et al. (1974) 
are the most acceptable values. 



Figure 1-2. The absorption index of ’esert 
aerosols in the solar visible region. 

K: Desert dust - Lindberg and Leude (1975) 

3K: 3X desert dust - Lindberg and Laude (1975) 

KK: Desert dust - Kondratiev et al. (1974) 

KJW: Average desert dust - Joseph and Wolfson (1975) 

KL: Limonite - Egan and Becker (1968), Kondratiev (1972) 

(2) The Solar Spectral Region: 1.1 pm <. X < 2.5 pm 

Lindberg and ^mith, (5974), and Kondratiev, (1972) presented the most 
recent data on the properties of desert aerosol clays in -this spectral region. 
Their data show a continuous extinction with an extinction coefficient between 
.01 cm ^ and 1 cm leading to 8 x 10 ^ £ k £ 2 x 10 ^ . Superimposed on this 
are band-absorption features centered on the water-vapor absorptiou bands. 



-4 -3 

The value o£ k in these bands may reach lx 10 < k < 1 x 10 and 

.002 < k < .02 in the 2.7 pm band. The absroption in these bands depends 

*\/ 'Ki 

on the state of hydration of the clay sample. 

We recommend the averaged values of n ■ 1.54 and k ■ .008. These values 
are based on spectral averaging using the solar spectrum as a weighting func- 
tion. The value of k may be twice as large (Brinkworth, 1972). 

(3) The Thermal I'.ifra-red Spectral Region: X 2. 2.5 pm 

All materials that may be present in desert aerosols have strong ab- 
sorption bands in this region and they may be bi-refringent. Peterson and 
Weinman, (1969), and Steyer et al., (1974) present results for quartz, and 
Lentz and Hoidale (1974), present results for silicate clays and carbonates. 

if . 

'■ ' ( 

The imaginary part^of the refractive index varies over several orders of 

• •" ; i 

magnitude in the shortwave and the 8 pm < X < 12 pm parts of this spectral 

I 

region. The maximum extinction coefficient may be 15-20 times larger than 

t 

in the solar range; Day et al., (1974). The same authors show that the 
determination of the imaginary part of the refractive index by transmissometry 
of samples dispersed in a matrix under-estimates the absorption at the peak 
by about a factor of three due to saturation. However, the spectrally inte- 
grated absorption as evaluated by Hie calculations based on the conventional 
data, gives the absorption within 15 percent. It follows thus that a "band 
model" for aercsol absorption based on these Mie calculations should have 
reasonable success in simulating the infra-red optical effect of aerosols. 

Evidence of the global homogeneity of the desert aerosol is demonstrated 
in Figure 1-3 which shows the imaginary part of the refractive index of desert 
aerosol as found in the Western U.S. by Lentz and Hoidale, (1974). Two similar 
sets of the spectral determinations of the same quantity from the desert of 
Israel are also shown. The large variability from dust storm to dust storm 
and even in one storm is also evidant from the difference in the two sets 
of measurements in the Israeli desert, which were taken by Fischer (1977). 
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Figure 1-3. Imaginary part of the refractive index of desert 
aerosols in the spectral region of the atmosphere's infra-red 
window . 

The real part of the refractive index may thus be taken as n • 2.45 
and the imaginary part may be taken from Lentz and Hoidale's model of the 
optical constants of the desert aerpsol. 

1-2-f) Size Distribution of Desert Aerosols 

The size distribution of desert aerosol has been investigated by Junge, 
(1972), Bullrich, et al. (1973), Lentz and Hoidale (1974), Volz (1970), Joseph 
and Manes (1971), Joseph et al., (1973), Joseph and Guetta (1975), Prospero, 
(1968), Carlson and Prospero, (1972), Peterson (1968), Bryson et al. (1964), 
Bryson and Peterson (1968), Lindberg and Laude (1975) Lentz and Hoidale 
(1974), Grams et al. (1972), Shaw et ai. (1973), Kondratiev (1972), Kondratiev 
et al. (1972), Schutz and Jaenecke, (1974). 
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On the average, the experimental 3ize distributions may be represented 
by a sum of two incomplete gamma-functions, one peaking between 1 and 10 ym, 

the other at radii less than 1 ym 
n(r) - n t (r) + n 2 (r) 

- a^r Ql expC-b^ 1 ) + a^ 2 exp(-b 2 r Yz ) 1-2.1) 

a 1 » 3.5 x 10 5 , cm 3 ym 3 = 2 . , b^ “ 20. ym **, * *50 

a 2 » 2.9 x 10 5 cm~ 3 ym 11- f a 2 « 10., b 2 » 10. ym J ', y 2 » 1.0 

Assuming the mass density to be 2.0, the corresponding mass is 
3 3 

1.25 yg/m and is 299. Ug/m for a "normal" desert aerosol. 

1 The n^(r) is introduced because the desert clays cannot absorb solar 
radiation in the amounts experimentally found. The simple way of providing 
for the absorption is to add a small mass of limonite in the form of snail 
particles. Kondratiev (1973) also found that limonite is only in aerosols 
with r<l ym and that the larger aerosols are clays. 

Not enough data are available to say much about the change in height 
of the size-distribution. Due to the efficient mixing by convection of the air 
over deserts we shall assume that the size distribution is independent of height. 

An estimate of the parameters of She size distribution of desert aerosols 
can be obtained in the following way: 

The refractive indices, as presented previously, have been used to 
compute the Mie optical depths for extinction in the solar spectrum, 0.4 ym 
< X < 1.1 ym, ard in the 8.0 ym < A < 12.5 ym thermal infra-red region for a 
range of values of the parameters of the size-distribution. In the infra-red, 
n 2 (r) is the main contributor to the optical depth whereas n^(r) and n 2 (r) are 
both important in the solar visible and near infra-red. 

Our comparisons show that the parameters cited for n 2 (r) u ield 
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a best fit between the optical depths computed from Hie theory in the 

8.0 pm < X < 12.5 pm region with those reported by Lentz and Hoidale (1974). 

We then take the optical depths in the visible, based on this best- 
fit size distribution, n^fr), and combine it with a set of extinction 
optical depths for llmonite with a size distribution n^(r). The parameters 
of the size-distribution n^(r) are then varied to obtain a best fit to the 
spectral extinction of solar radiation. 

The best fits, n^(r) and n 2 (r), are compared with various size distri- 

i ! 

butions in Figure 1-4. The sum of these two distributions constitute our 

i i 

I I 

proposal for the avor^ge size distribution of the desert aerosol as found 

/ 1 

1 ! 

over the oceans or over deserts during calm conditions. 

I 

1-2-g) The Maximal Amounts of Desert Aerosol 

In order to describe the extreme effects of desert aerosols on thermal 
infra-red heating rates as measured by Bryson et al., (1964), for example, 
quite large optical depths must be assumed (t ^ 3 - 6) . Above "normal" optical 
depths are also needed to partially explain the decreased cooling in the lower 
troposphere observed by Peterson (1968) . This means that between 10 - 20 x 
normal amounts of 112 (f) are needed to explain the Bryson et al. case, and 
1 10 x normal are required to account for the Peterson measurements. Measute- 

ments of the concentration of large particles over the Sahara Desert by SchUtz 
and Jaenicke (1974) indeed show concentrations of large particles which are 
5-10 times greater than were assumed in the derivation of our "normal" desert 
aerosol; 

The average optical depth at X = 0.5 pm of the aerosol during Khansinic 
phenomena over Israel is .4, however, optical depths exceeding 1 occur about 10 
percent of the time, Joseph and Manes (1971), Joseph, Manes and Ashbell (1973), 
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Joseph and Wolfson (1975), Grassl (1974). The data on visibility over the 
Middle East collected by Lentz and Hoidale (1974) show that visibilities 
less than 11 km occur about 50-150 days every year. This implies that the 
visible horizontal extinction coefficient is larger than .35 km If 
typical aerosol scale heights lie between 1 and 3 kms, a vertical optica) 
depth between .35 and 1 occurs about If'-' 1 days every year. Carlson et al. 
(1973) report a reduction of global solar radiation equal to that of thin 
clouds and a Linke turbidity factor of .3 at X = .0.5 pm during the occurrence 
of dust over the Atlantic. This value is equivalent to a vertical optical 
depcb of .69. 

The range of masses of the desert aerosol modelled in this study is 

3 

therefore assumed to vary between 30 and 3000 Pg/m at the surface. We shall 
denote a normal desert aerosol by BF = 1, signifying a surface mass-concentra- 
tion of 300 pg/m 3 . Aerosol scaling factors of BF of 0.1 and 10 denote "light" 
and "heavy" desert aerosols respectively. 


1-3 Models of the Clear Atmosphere 


Two models were assumed for the clear atmosphere. The first is that of 




a tropical atmosphere defined by McClatchey et al (1970), and the second 
is that given for April 24, 1966 over New Delhi, India by Peterson (1968). 

The surface temperature was assumed in each case to be equa.' to that of the air 
at the surface. The main difference between the two is the higher air tempera- 
tures below 700 mb in the New Delhi case. The air temperatures at the surface, 
at 900 mb and 800 mb, are about 10°K higher. That at 700 mb is two degrees 
higher. The water-content is about the same. The ozone-profile from the 
McClatchey model was used in both cases. 

The optical properties of water-vapor at 8 pm < X < 13 pm were taken 
from the study of Cox (1973) and both continuous and e-type absorption 
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Figure 1-4 Size distributions of desert aerosols (a) Rajasthan desert - 
Peterson (1968), (b) Colorado - Whitby (1974), (c) Texas - Gillette et al. 
(1972), (d) Rajasthan desert - Bryson et al. (1964), (e) Sahara sand 
storm - Schlitz and Jaenicke (1974), (f) Camp Derji Sahara, sand storm - 
Schlltz and Jaenicke (1974), (g) Mean of Sebha, Sahara - Schlitz and 
Jaenicke 0.974), (h) Atlantic Ocean - Schlitz and Jaenicke (1974) f 
(i) Grams et al. (1974), (J) and (k) de. uisi, private communication. 
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were taken into account. All calculations were done in spectral intervals 
of AX a 0.5 ym in that spectral region. The data used were linearly inter- 
polated from the data summarized in Table 1-2. Temperature and pressure 
dependence of the absorption coefficients were taken into account as shown 
in Equations 1-3,1 to 1-3,6. 


T.blel-1 I Optical Properties of Ceseous Ateosph.rlc Absorber' 


(a) Water-vapor 


Vevelength, (vie) 8.00 8.15 8.70 9.10-9.55 10.00 10.55 11.10 11.80 12.50 11.35 

K (293*K), j/cb 1 # 090 .080 .074 .069 .062 . 068 . 090 .100 .302 .108 *113 

** (293*10 lg/c m*f 19.80 17.00 14.50 12.00 9.30 9.30* 9.50 *10.00 12.00 14.00 20.00 

» 

/{ , (b) Ozon a 
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0.000 
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0.000 0.000 

0.000 
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0.000 

.008 .020 

.051 


The ozone amounte are calculated, in units of atm-cm, reduced to 


S.T.P. 




1-3,1 


(McClatchey et al 1970) . 

The numerical factor reduces the ozone amount to S.T.P. 
£ = p/p o 

-2 

p is surface pressure in dyne cm 

-l 

g the gravitational constant in cm sec 
c the ozone concentration in g/m^ 
p the air density in g/m" 
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Water-vapor amounts were modelled by 

tZ 
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K_0, e-type “ .622g 
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1/2 

l f (A.) dJL 
P(I) 


S, f^f9\ Al 

C(293 # K) 1 W d * 


1-3,2 


1-3,3 


6 —2 

where m 1.01 x 10 dyne- cm 

q(l) - 1.63xl0~ 2 /g 
-2 

g - 981 cm sec 
T x - 273°K 
T ■ temperature, °K 


2 1/2 

C(T) “ 2 + x(T) - [(2 + x (T)) 4 - 4) 4 ' 4 


1-3-4 


x(T) 


■622 exp { - Cl - (c 2 /T) } 
q(£) f(A) 


1-3-5 


[ £ 3 * 56 

f(JL) (.130) 

[ ° 

13.01 


3.56 


.130 < l < 1 
.001 < Z < .130 
l <.001 


c 2 - 2878. 2’K 


(Cox 1973) 


1-3-6 


The models for absorption by gases outside the 8.0 pm < X < 12.5 Mm 
window are based on those of Rodgers and Walshaw (1966), and Rodgers (1967) 

and used conventionally in the manner cited by R. Bursztyn (1974), to evaluate 

fluxes in 65 pressure intervals from the surface to 1 mb. 
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2-1 He ating Rates due to Desert Aerosols 


2-1-a) General Discussion 

The role of aerosols in the general circulation of the atmosphere and 
in climate has been the subject of extensive discussion. The problem of 
aerosols arises in two ways: 

The first problem is the determination of the impact of aerosols on 
atmospheric and surface heating. The second problem is the determination 
of their indirect role on the energy balance of the atmosphere through their 
effect on the formation and the properties of clouds. We shall concern 
ourselves solely with the first problem. 

2-1-b) Heating Rates for Numerical Models of the General Circulation 

Random errors in heating-rates have little effect on General Circulation 
Models according to Joseph (1966), and Washington (1971). It is important, 
however, that the heating rate departs systematically from a suitable stan- 
dard and that its trend is known. Systematic changes in radiative properties 
of the atmosphere have quite significant effects, see Joseph ibid, Fels and 
Kaplan (1974), Manabe and Wetherald (1975). For example, systematic pertur- 
bations of the heating rate greater than 0.5°K day * affect the NCAR G.C.M. 
according to Washington (1974). 

An important requirement that any radiative heating rate scheme must 
satisfy to be useful in numerical modelling of the atmosphere is speed of 
computation. We address ourselves now to the development of fast and rela- 
tively accurate radiative transfer model to compute the heating rate of the. 
desert aerosol in the atmosphere. 



2-1-c) Heating Rates in the Solar Spectrum due to Aeroa 


Several studies of solar heating due to aerosols are available 1 i the 
literature e.g. Yamamoto et al. , <1974), Dave and Braslau, (1974) and Braslau 
and Dave (1973, (a), (b)). 

An important problem in modelling heating by an aerosol is the variation 
of its mixing-ratio with height la the atmosphere. This leads to the need 
for many layers in order to adequately model the variation of the albedo fqr 
single scattering with height. 

A second problem is the spectral integration of fluxes and heating 
rates. The extinction of the aerosol usually decreases slowly with wavelength 
and it is small. In the appropriate spectral ranges, compared to the absorp- 
tion b> water-vapor and other gases' absorption bands. Rayleigh scattering 
is small compared to absorption in the solar infra-red. We therefore divide 
the solar spectrum into two regions. One contains the range .3 pm <A< .8 pm 
and all "windows" between gaseous absorption bands; aerosols are assumed to 
be optically active and Rayleigh scattering by molecules must also be con- 
sidered. The second region is that of all gaseous absorption bands; there 
we assume that only gases are optically active. 

2-2 The 6- Ed ding ton Approximation 

A simple model for atmospheric radiative heating rates due to the 
presence of aerosols will now be developed : 

The direct use of the Eddington approximation in calculating solar 
fluxes and heating rates leads to large errors for optically thin aerosol laden 

atmospheres. The reason for this is chat single-scattering is usually dominant 

* 

due to the relatively small total optical depths, T < 1. In this case, 
the actual phase-function is strongly peaked, while the Eddington approxi- 
mation assumes that the peak in the scattering function is relatively small. 

The Eddington phase function is 

1 ♦ (“jcose, 2-2,1 
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The forward Co backward ratio Is thus 
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r = -L(Q!I. ,13 

“ F(180 # ) 1 - « 


2 - 2,2 


Van de Hulst (1971) showed that <o^ • 3g, where g ■ P(0)cos0 d(cos0) 

-1 

is the as;itmetry factor. In the solar spectral range, g ~ .7 for aerosols. 
Therefore, R ~ - 2.8 while its actual value may be of the order of a 1000 
and positive . 

This difficulty can be reduced if the forward peak of the phase 

function 1 b replaced by a suitably weighted Dirac delta-function. The phase 

! • 

function can be approximated by: 


P(M) - 4irf<S(y - 1)6{$) + (1 - f) (1 + 3g’u) 
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where y “ cos0, f is a weighting factor defined belcw and g* is the 
asymmetry-factor of the truncated phase function. 

The application of this approximation of the phase-function leads to 
a scaling of the optical depth, albedo for single scattering, and the 
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The factor f may be found by 
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Then Che second moment: 


PjCcosejdn = f 


2-2, 5c) 


where P^Ccosd) is the second order Legendre polynomial. 

Any phase-function may be approximated by a Henyey-Greenstein function. 
Van de Hulst (1971) showed that the Henyey-Greenstein function may be expressed 


P(0) - l (2£ + l)g £ P.(cos0). 
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Therefore, it follows that 


f - 8 
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We can derive fluxes from an Eddington approximation which employs the 

transformed quantities defined in Eqs„. 2-2,4 and 2-2.7. The quantities I 

' o 

and 1^ are given by Shettle and Weinman, (1970). 


X - C. exp(-Kt’) + C exp(Kt') - a exp(-x'/w ) 
OX c o 
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In order to teat the accuracy of this method, the fluxes and flux 
differences were extensively compared to similar calculations using doubling 
method, see Hansen and Travis (1974). Table 2-1 shows sample comparisons 
for the moat stringent demands on the 6-Eddington method, namely small 
optical depths, t* * .0282 and .282, a surface albedo of zero and a range 
of solar zenith angles 1.0 > > 0.2. All fluxes and flux differences are 

given in units of the normal incidence solar incident flux at the top of the 
aerosol layer, T “ 0. 

The total downward flux, computed using the 6-Eddington approximation, 
is usually equal to that evaluated by the doubling method to two or three 
significant digits. The same can be said for the diffuse upward and the net 
fluxes. The net flux differences, which are proportional to the heating- 
rates at any level, indicate a systematic error in the 6-Eddington method. 



































































For the smallest optical depth considered, .0282, it is -15 percent and 
-11 percent for the upper and lower half of the aerosol layer respectively. 
In the case of an optical depth ten times larger, namely .282, the same 
errors are -22 percent and -7 percent respectively. The largest error among 
those enumerated above, namely -22 percent in the upper half of the aerosol 
layer with a total optical depth of .282, leads to an error in heating rate 
of -.07°K dcy" 1 . 

It is thus clear that the 6-Eddington method presented in thio study 
is adequate for the calculation of monochromatic fluxes and heating rates 
In an absorbing and scattering layer of aerosol. 

We now treat the parameterization of the monochromatic h‘ a ting in an 
atmosphere containing absorbing and scattering aerosols embedded In an 
ambient Rayleigh-scattering atmosphere. 

The scattering phase-function for the mixture of a Rayleigh-scattering 
molecules with a Mie-scattering aerosol, may be expressed by n weighted Bum 
of the separate phase functions, see Deirmendjlan (1969). 





is the scattering coefficient of the aerosol; 
is the scattering coefficient of the molecules; 



and: P(9), P(G) are the phase functions for the aerosol and molecules 

a m 

respectively. 


The asymmetry factor of the mixture is: 
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Because: 


g„ 


= 0 
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Consequently, the amount of energy scattered into the forward peak is: 
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The albedo for single scattering of the composite medium is: 


ui (1 + x(h)) 
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, Where u is the albedo for single scattering of the aerosol. The 
a 

molecular scattering coefficient s ^ is proportional to pressure: 


(h) „ s (o) h / H } 
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| The height profile of the aerosol extinction roefficiert may be approximated 

I 

jj by an exponential profile with a scale height ■ 1.0 tan for the case con- 
| eidered by Dave (1975). 
i a » 8 ^ exp {-h/H„} 
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We now parameterize this composite atmosphere by evaluating a pressure- 
weighted average albedo for single-scattering u and a symae try-factor g : 
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These equations, were evaluated in the two spectral regions 0.3 ym <X< 0.8 y% 

/! ' 1 

and 0.8 ym <X< 2.5 ym. Only in the first region is this correction of the 


scattering parameters of importance. 


Table 2-2 summarizes the results oi this approximation for the scattering 

i 

parameters for <• mparison to Dave and Braslau's (1974) work. The molecular 
scattering coefficient at the surface was taken to be .019 km 1 


Table 2-2 The Effect of Molecular Scatte*ing on an 
Aerosol Laden Atmosphere in Spectral Range 
0.3 ym < X < 0.8 ym. 


Aerosol Only 

Aerosol + Molecule 

Dave Cl (1974) 


8 (o) - .096 km' 1 

a 

km" 1 

a m 

to ** .910 

a 

w - .978 

g a - .-oo 

g - .143 

f - .490 

a 

7 - .021 
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The results of this comparison are shown in Table 2-3. 


tab la 2-3 “ Coaparlooo of C- Eddington ' Computation* with Tboae of Daea (1974) 


Solar Angle 


e 

o 

degree*) 


7.5 (Binds) 


Ket Flux {Net Flux {Net Flux {N^l flux (Wet Flux 
Bottom 


Total Absorption Totel Absorption (Dove Cl 
l 1 



.8 - 2.5 pm 


•8 - 2 . 5(bxnds) 



Flux Up Flux Dovn Flux Up Flux Down Flux Up Flux Dovn Rayleigh Present Cave (Cl) 
Top Top Top Top Top Top X X X 

(v-o~ 2 ) (v-o~ 2 ) (v-e 2 ) (v-* 2 ) (w-n 2 ) (v-n 2 ) 


<«5 Reflection* Surface Albedc 

- 0 

82.6 

799.7 

.278 

64.72 

82.4 

772.5 

.294 

62.52 

81.8 

692.6 

.336 

56.06 

76.0 

399.9 

.483 

32.36 

53.5 

136.9 

.776 

11.24 


15 Reflection. Surf see Albedo 



17.3 

1 

» 64.7 

1.9 

!n.2 



Reflection. Surface Albedo 


207.7 

799.7 

147.6 

138.9 



527.6 

2) 

91.6 

33 
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Table 2-3(a) shows the results for three separate spectral 

ranges .3 Mm <X< .8 Mm, .8 Mm <X< 2.5 Mm between the gaseous absorption bands 

and .8 Mm <X< 2.5 Mm in the gaseous absorption bands. Spectrally averaged net 

2 

fluxes are expressed in units of watt/m , at the top of cur atmosphere and at 
its bottom. From this set of data we calculate the total absorption in units 
of the incident solar flux to be compared with the exact results of Dave. The 
gaseous abosrption was modelled by fitting a "best-fit" curve to the absorptivity 
given by Dave (his Model B) , as a function of the solar zenith angle. Ve used 
the water-vapor profile given by Dave. The resulting equation for the combined 
transmissivity of all absorbing gases in a Rayleigh scattering atmosphere 
for water-vepor amounts up to 8 g/cm and for the whole solar spectrum, is 


where 


. * 1/2 

T(~) * .94 exp { - .063 (^-) } 

"o W o 


2 - 2,21 


f 

uMt^.tu) « I 


p < h > (J^)dh. 

water p(o) 


is the "reduced" water-vapor mass between two altitudes when the sun is 
at a zenitn angle *->s The transmitted intensity at a given height 


F< f> ‘ Vo T( f> • 

o o 


2 - 2,22 


where the incident solar flux is that over the whole spectrum. This trans- 
mission function is similar to many other empirical equations derived in the 
literature, see e.g. Kondratiev (1972) except for being a transmissivity for 
the whole solar spectrum and the inclusion of the effect of Rayleigh scattering 




- / 
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on the absorption. It is interesting that a square-root dependence ia found, 
even though we modelled the total absorption by all gases throughout the solar 
spectrum by Just that of the water-vapor. 

The last step in the approximation of absorption by gases is the 
assumption 


0**‘)V o a - ) - I 0 U 0 U - T(^) ) 2-2,23 

I U O . 

I 

or that the total energy, throughout the solar spectrum, absorbed by the 

I 1 ; f 

gases may be thought of as being concentrated in the absorbing bands of 
water-vapor at .8 yu £ X £ 2.5 ym. 

: \ 

The net fluxes in the water-vapor absorption bands are calculated 
with- the help of Eq. 2-2,17 by using the following equations 


A ^ 

F(top)- - I o (bands )y o [l - T(J-)T(1.67u*)A] 
bands *o 


2- 2, 24a) 


F (bottom) “ - I (bands) y T(~~) fl - A] 
, , o o y 

band 8 o 


2-2 ,24b) 


where A is the surface albedo and I o (bands) - .4255 I Q . 

On comparing the last two columns in Table 2-3, (a), (b) end (c) , it is 
seen that atmospheric absorption is accurate enough to be used in a general 
circulation model. 

General circulation models also need fluxes at the top and bottom of 
the atmosphere to determine energy balances. We therefore also compare 
the reflection by the aerosol-molecular atmosphere with that of the 



i 


5 

} 

i 

is 

7 


e 

i 




more exact results. The results are shown in parts (d) , (e) and (f) of Table 2-3. 


The application of this method yiel/3 albedos that are in good agreement 
with those obtained from more accurate models as is shown in the last two 
columns of Table 2-3 (d) , (e) and (f). 

The same ideas may be applied to the calculation of the downward flux 
at the surface. A comparison of the 6-Eddington approximation to the "global" 
flux at the surface with the "accurate" results is given in Table 2-4. 

A 

Table 2-4. The Downward Fluxes at the Surface in Units of (J I " Comparison 


with Dave and Braslau (1974) and Dave (1975). i 



The comparisons of our computations with mere exact computation shown 
in Tables 2-2, 3 aud 4 indicate that it is possible to describe the effect of 
aerosol on the solar heating rate of a clear atmosphere in a way that is 
simple, fast and with a predictable systematic error not larger then a few 
percent. 


fUL 
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We now address ourselves to other physical errors In computing heating 
rates in an atmosphere containing a mixture of absorbers and acatterera 
throughout the solar spectrum. We defined three spectral intervals 
.3 pm < X < -,8pm, .8 pm < X < 2.5 pm between bands, and .8 pm < X < 2.5 
in bands, and averaged the effect of the Rayleigh atmosphere. 

The division into two main spectral regions does not result in large 
errors, see Joseph (1966), (197ia) ,Lacis and Hanson (1974), and Dave (1975). 
The neglect of ozone in the troposphere is not serious. The vertical optical 
depth of ozone at X^O.^ pm, the center of the Chappuis band is of the order of 

I .04 for normal ozone amounts (.384 atm-cm) and more than 90 percent of the 

j ,! t 

ozone is above 1C km. | The effect of ozone in the stratosphere on tropospheric 

i • I' i 

j fluxes may be modelled (by decreasing the solar flux at the tropopause by the 

! ' ! 

ozone transmission, see Joseph (1966). 

Molecular scattering in the spectral region 0.8 pm < X < 2.5 pm can be 
neglected because it i6 insignificant. The fact that we do not account for 
absorbers other than water-vapor in this latter spectral region does not lead 
to large errors, e.g. Figure 8 of Braslau and Dave (1973) and Table 2-3. 

The neglect of the over-lap of the water-vapor and aerosol spectral 
characteristics in the spectral region .8 pm < X < 2.5 pm will lead to errors 
of about 2X in the heating-rate. The error is due to the over-lap in the 
wings of t'le water-vapor bands and is small because there is little solar 
energy in this optical range. 

The heating rate is evaluated directly through the use of 


1 -T’/P 0 

g? - And - u>’) U o (t’) + ^e ] 


2-2,25 


3t (h> V A; " “ p(n) 




2-2,26 




The calculation* of the heating-rate have been made by using an 

A 

optical depth, weighted by the solar apectrum, l o (X) , 

. 8ym . • 6pm 

x ■ f t(X) I Q (X)dy J I Q (X)dX 

.3 Vim ' .3pn 

in the range 0.3 um < X < .8 vim, and 


"average" 


2-2 »27a) 


t » E 




U, ttt)I 0 (X)dX 


A k 


I (X)dX 
o 


2-2 ,27 b) 


where the AX are all spectral intervals in the solar spectrum between .8 and 
2.5 vim where there is no gaseous absorption. (.80-. 88, 1.00-1.10, 1.18-1.28, 
1.52-1.68 pm). 


Models of the general circulation are insensitive to errors in the 
heating-rates wh^-ch are less than .5°K -day \ Washington (1974). The 
error in our approach is acceptable whereas the additional time required 
to rigorously compute the heating-rate in an atmosphere with two mixed 
scatterers would not nermit the inclusion of aerosols in a model of the 
general circulation 


2-3 The Atmospheric Heating lute Models in the Solar Spectral Range for 
Application to G.C.M.'s and Climate Models . 

The technique outlined in the previous section has been applied to 
models of the desert aerosol over land and over sea. 

The optical characteristics at the surface are ohown for a "normal'' 
desert aercsol in Table 2-5 es used in our model. 


f 


Table 2-5* Optical Characteristics of "NORMAL" Desert Aerosols 

A a A A 

I q • 1392 w/m , Solar Constant ( 61 ^/ 1 ^ in percent. 


AX [pm] 


.3-. 4 .4-. 5 .5-. 6 .6-. 7 .7-. 8 


8 (0) [km A ] .230 .210 .196 .175 .155 
a 

a) .934 .958 .970 .975 .976 
g .670 .670 .670 .670 .670 
61/1 (I) 7.8 14.4 13.8 11.8 9.6 


.3-. 8 .8-2.5 

.193 .055 

.970 .920 

.670 .670 

57.4 42.6 


Several model distributions with height of the desert aerosol were 
used and these are shown in Table 2-6. 

Table 2-6: Desert Aerosol Distributions as a Function Height 


Top 

[km] 

Scale Height 
[kr/J 

Notes 

7.3 

1.31 

Overland - Normal 

7.3 

2.8 

Overland - Dust Storm 

5 

OO 

Over Ocean;* 8 (h) « ( 

a 


*This distribution is required as input for G.C.M.s; however, lidar 
parameters are computed from the overland normal model only. 


The inclusion of the effect of Rayleigh scattering on the heating rates 
over land is similar to that applied to the aerosol model Cl of Dave. 

The results are shown in Table 2-7 (a) for 0.3 pm < X < 0.8 pm 
and should be compared to those in Table 2-5. 


i \ 




Table 2-7: Optical Characteristics of Atmosphere Containing Desert Aerosol 


(a) Overland Light (BF » .1) 

Normal (BF «* 1.) 

Heavy (BF 

p (o) + gCo)^ -1 ) .038 

a m 

.211 

1.949 

Z, .987 

.974 

.970 

g .176 

.491 

.645 


(b) Over sea 
(In layer) 


g(°) + g(°)(k m ~ 1 ) 

a , m 

.031 

.205 

1.942 

5 j 

-.981 

.972 

.970 

5 ii 

’.098 

l 

.153 

.162 


[ 

Note that BF is defined In section l-2g) . 


The analysis of the desert aerosol aloft, in a layer between 

i 

2 and 5 km, involves the same averaging technique as before, but applied 
only in the layer itself. The results are shown in Table 2- 7 (b) . 

The more difficult problem is the adaptation of the boundary conditions, in 
the 6-Eddington solution necessary to take into account the existence of 
Rayleigh scattering above and below the aerosol layer. This is handled in 
the following way: 

The downward flux of solar radiation at the top of the aerosol layer at 


5 km is modelled by 

„ -T (5)/u -1.66T (5) 

, , . _ . m o,l — e m , 

Fi(5) = I o h o [e + 2 1 


2-3. 


~ -.07/p 

“ Vo le ° + * 051 

where the solar flux incident on the atmosphere is decreased by the direct trans 

-t_(5)/y . . 1 - e" 1 * 66T m (5) 

mission, e .and enhanced by the diffuse transmission, 2 * 

through a Rayleigh atmosphere down to 5 km. 
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The lower boundary condition, now applied at 2 km, has to take into 
account the existence of the Rayleigh scattering between 0 and 2 Inn. The 


upward flux at 2 km is approximated by 


Ff<2) s [M o I o ex P (-T’(2,5)/u o > {exp( - 


AMO, 2) 


1 - e 


-1.66 At (0.21 
m 


+ tt [I (2) + 4 M2)] 

O 3 1 


-1.66AT (0,2) 
m 


rl + e 


-1.66 Ax (0,2) 

m 


where 1^, I q ( 2) arvd 1^(2} include the corrected upper boundary condition 

(Equation 2-3.1) and At m (0,2) is the spectrally averaged vertical optical depth 

of the Rayleigh atmosphere between the surface and 2 km, namely .03, A is the 

surface albedo. An analysis of the first two expressions in curly brackets, 

involving At , shows that one may equate the two with an error of less than 
m 


.03, and express the l-'wer boundary condition by 

Ft(2) ° F+(2) (.975) 2 A 

Effective $ <- 575>2i - 
where F*(2) is the downward flux at 2 km. 


2-3.3 


2-4 Results - Desert Aerosol Over Land 


Atmospheric solar heating rates due to desert aerosols were evaluated 
as a function of height for various values of the solar zenith angle, and 
the surface albedo in each of the two spectral ranges (.3 pm < A < .8 pm, 

.8 pm < X < 2.5 pm) . The height increments were 0.5 km and 1.0 km. The 
heating due t« wnter-vapor in the same tropical atmosphere was also available, 
as was the total heating. 

We shall only present a few examples of our results: 

Figure 2-1 compares the ^eating rates due to aerosol in the solar spectrum with 
that due to water-vspor . The atmospheric profiles of temperature and humidity 




Figure 2-1. Comparison of solar heating rate due to aerosols with that '' 

due to absorption by water vapor in the interval 0.3 ym < X < 2.5 ym. ‘ 

Water vapor absorption function is from Kondratiev (1972). Numbers ~j 

identify the total optical thickness, T*. Surface albedo, A ■ .15, j 

•J y - .6, H - 1.3 km. | 

° a •{ 

5 

' ' are taken from Peterson (1968) and they are typical for a desert atmosphere '■ 

• ? 

! over India. Increasing total amounts of desert aerosols are indicated by . ' 

the optical depth, T*, at the surface. The curve marked x* “ .395 is for an ; 

! I 

i 

average amount of desert aerosol in a clear atmosphere. The average amount ; 

i ' i | 

of desert aerosol leads to a heating tnat is about one half of that due to 

j*“ • 

! water-vapor at the surface and decreases with height due to the aerosol 

scale height, namely H = 1.3 km, vhich is for the desert aerosol over the 

a 

. i 

Indian desert. A dense desert aerosol may have a vertical optical depth 

• j 

larger than one. The curve marked 1* = 3.95 illustrates probably the ex- , 

treme effect of a desert aerosol. In the lower troposphere, the heating 
due to this amount of aerosol is much larger than that due to water-vapor 
and a peak in the heating rate develops aloft. 
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P.-5 Desert Aerosol Layer Aloft In a Moist Tropical Atmosphere over the Ocean . 

A layer of desert aerosol between about 2 and 5 km over the tropical 
Atlantic Ocean may significantly affect the transfer of solar radiation, the 
heating-rate profile, and the stability of the atmosphere. We have therefore 
modelled the effect of a desert aerosol layer aloft over the ocean. As an 
example of our results we show here in Figure 2-2 the results for a thin 

(BF = .1) and a normal (BF ™ 1.0) desert aerosol. The optically thin model 

I ; 

aerosol, with an optical depth at X <* .5 pm of 0.06 has a constant mass- 

, i 

i ' '3 

concentration of 30 pg/ra in the j.ayer. For a normal desert aerosol, these 

i ! 1 

I f 

quantities should be multiplied by a factor of 10. The ambient tropical 

.. /; ■ [ 2 

atmosphere has a water content of 3.57 g/cm and a mixing-ratio at the sur- 

i ! 

face of ,016 g/g. 

It is seen that a thin desert aerosol layer aloft has only a minimal 
effect on the solar heating rate. A normal (BF - 1) desert aerosol Increases 
the heating in the layer by more than .5°K day ^ at its bottom aid by 
• 8? J K day * at its top. This is a very large change in heating, which will 
strongly affect atmospheric stability. 

In Figure 2-2 Part (b)we compare our theoretical Increase in solar heating 
due to the aerosol over an ocean with the difference between measured and 
calculated clear-air 6olar heating rates, as given by Reynolds et al (1973). 
Unfortunately, the measured data have a large scatter, and comparison of 
individual days is difficult. The comparison with our results shows that 
in the aerosol layer, the theoretical differences between the effect of 
layers containing 30, and a 300 pg/m'* cover the measured differences 

3 

(90 to 160 pg/u ). It is seen that the experimental differences reach a 
maximum near 3 km where the presence of aerosol is usually detected. The 
difference below 2 km is near zero and in the lowest kilometers chere is 



another relative maximum, which might be due to marine haze or sea-spray. 


38 



Figure 2-2. Solar heating profile for desert aerosol aloft over the ocean. 

3 

(— * — *-) Computed BF *= 0.1; T (2) * .06, mass concentration 30 pg/m BF = 1.0, 

3 

x(2) = .58, mass concentration 300 pg/m ( ). Tropical atmosphere 

[McClatchey et al. (1971)], A = 0.05, P q “ 0.6. Part b: Difference between 
heating ir. Dusty and Clear Atmospheres. Mean of eight observed profiles, 

KMS eiror i 0.5°K day ^ ( ). 

it thus seems that a desert aerosol aloft over the ocean may signifi- 
cantly affect the distribution of solar heating. 

2-6 The To tal Daily Solar Radiative Heating Over a Desert 

The total daily solar heating over a desert is a quantity of interest 
for a model of desert climate or of the role of the desert in the general 
circulation of the atmosphere. 

We modelled this quantity by integrating the solar heating of 
a typical desert over the hour angle of the aun for the case of a desert 
aerosol with a scale-height of 2.8 km distributed up to 400 mb in a clear 
atmosphere. The surface albedo was A ** .25. 

The results are shown in Table 2-8 for the dust, the water-vapor and 
for their sum for three different amounts (BF ■ .1, 1.0, 13.0). 
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For a thin layer (BF - .1), the heating due to water is dominant. A 

S I 
‘ ( 

| normal desert aerosol (BF - 1.0) will add approximately 50 percent of the 

heating due to water up to a height of three kilometers and be slightly 

less significant at higher levels. The very dense aerosol case described by 

I 

■ i BF - 10. will more than triple solar heating. 

' I 

I 

( 

i 

i 


{. — • — ■ 

ht t 
(km)\ 

Table 2-8 . The Daily Solar 

H o 0 > Aerosol Heating Rate 

teatindBF - . 1 BF - 1.0 BF « 10. 
Rate/ 

H 

eating 

TOTAL 

Rate 

(b:- 

Over 

1) 

a Desert 
TOTAL (BF- 

(°K day -1 ) 

1.0) TOTAL (Bf-10.) 

0 

.63 

.04 

.33 




.67 


.96 


2.28 

! * 5 

.60 

.03 

.28 




.63 


.88 


1.82 

; l.o 

.57 

.03 

.26 




.60 


.83 


2.40 

| 1.5 

.54 

.02 

.23 




.56 


.77 


2.22 

2.0 

i; 

.52 

.02 

.20 




.54 


.72 


2.13 

»i ^ 

i! ^ 

.49 

.02 

.18 




.51 


.67 


2.01 

3.0 

.47 

.02 

.16 




.49 


.61 


1.89 

' 3.3 

i 

.45 

.01 

.14 

1.30 



. 46 


.59 


1.75 

! 4.0 
1 

.43 

.01 

. 13 




.44 


.50 


1.63 

| 4.5 

.41 

.01 

.11 

1.09 


• 

.42 


.52 


1.50 

! 5.C 

.39 

.01 

.10 

.58 



.40 


.49 


1.37 

'5.5 

.37 


.09 

.88 



.37 


.46 


1.25 

I 0.0 

.35 


.08 

.79 

! 

.35 


.43 


1.14 

: b - b 

. 3 4 


.07 

.70 

i 

.34 


.41 


1.04 

! 

..32 


. u 6 

. 63 

l 

i 

.32 


. 38 


.95 

I 7 - 5 

.30 


.06 

.57 



.30 


.36 


.87 

j 8.0 

.27 


.05 

.50 



.27 


.32 


.77 

j 8.5 

.22 


.04 

.45 i 


.22 


.26 


.67 

| 9.0 

.14 


.04 

.40 | 


. 14 


.18 


.54 

2 


! 

L 







1 


\ 


i 


40 


.. ' i 


2-7 Heating Rates In the Thermal Infra-red Spectrum 


2-7-a) Introduction 


i ; 


The calculation of heating rates in an absorbing, scattering and 
emitting atmosphere in the thermal infra-red poses some problems to the 
modeller. The main problem is that spectral fluxes and heating rates 

must be integrated over spectral regions in which the albedo for single 

i \ 

scattering varies markedly, see Sargent and Beckman (1973) and Uar.g and Domoto 

i! I 

(1974). |! Furthermore, the sources of radiation are both inside and outside of 

■ ' ! f 

the medium. An additional consideration is that the treatment of the fluxes 

( 

and heating-rate should be as close as possible to that accorded to absorbing 
gases, in order to avail one-self of the techniques developed over the years 
for treating the transmission of black-body radiation through the atmosphere. 
Our aim i'. this section is to develop a method of calculation that is adequate 
for use in a numerical model of the general circulation. 

2-7-b) The Radiative Transfer Equation in the Thermal Infra-red Spectrum 
The complete monochromatic radiative transfer equation in the thermal 
infra-red is 

^ (t, y, <J>) - I (t, y, $) - (1 - u>(c))b(t) 


^ : 


, 


<»(T) 

417 


P(y, <fr;y\ ♦’):-.(t, y\ ♦ , )dn' 


2-7,1 


471 


where all parameters have the usual meaning ai d B is the local Planck source 
function. 

The phase-f unction for aerosols in the thermal infra-red is somewhat 

peaked with an asymmetry factor .2 < g < .6, see Peterson & Wei-man (1969), 

A 

Deirmendj ian (1969), and Sergeant & Beckman (1973). We shall again reduce the 


* i 


/ 


/ 
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problem to a more isotropic scattering case by removing the forward peak from 
the phase function in the manner described in section 2-2. 

The radiative transfer equation is thus reduced to: 


"ft-* 


' ~ f 

4ir J 
4ir • 


P'(0)I dO’ - (1 - u’)B 


2-7,2 
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where the albedo for single scattering and optical depth are redefined by 
Eqs. 2-2,4, and 2-2,7. 

Table 2-9 shows typical values of the parameters w' and x* for three amounts 
of aerosol at X - 8 um and 11 pm. 




Table 2- 9. Typical Values of the Reduced Optical Depth and Albedo for Single 
Scattering in a Humid Tropical Atmosphere 
("top" of aerosol 400 km, scale height: 2.8 km) 


\ 

E 

Pressure 

Light 

Light 

Normal 

Normal 

Heavy 

Heavy 

1 


(mb) 

X* 

8ym 

A- 

llUm 

X- 

8ym 

X* 

llym 

A» 

8ym 


•lUm 


? 

T r 

(S' 

T* 

0)' 

T' 

w' " 

T* 

U)' 

T 1 

u» 

T ' 

w* 

i 

i 

1 

V 

1000 

1.25 

.004 

.78 

.004 

1.47 

.040 

.?0 

.033 

3.13 

.273 

1.92 

.216 

X . 
} 

1 

» 

r 

920 

.63 

.005 

.41 

.005 

.80 

.051 

.50 

.042 

2.04 

.323 

1.27 

.250 


\ 

840 

.31 

.009 

.21 

.0070 

.43 

.081 

.28 

.063 

1.27 

.42 

.83 

.311 

<' 

l 

4 

760 

.15 

.014 

.11 

.010 

.23 

.122 

.15 

.089 

.80 

.509 

.54 

.366 

) i 

{ 

[ 

680 

.07 

.024 

.05 

.016 

.12 

.189 

.08 

.130 

.49 

.597 

.34 

.420 

L* 

* 

r 

\ 

600 

.03 

.042 

.02 

.026 

.06 

.282 

.04 

.184 

.28 

.666 

.20 

.464 

I 

: 

520 

.01 

.074 

.01 

.043 

.02 

.400 

.02 

.253 

.14 

.716 

.10 

.498 

j 

- 

440 



.004 

.132 

.004 

.072 

.01 

.525 

.060 

.334 

.04 

.748 

.03 

.523 

* 

* 


Representative values for g' in the truncated phase function of the 
desert aerosol in the 8 ym < X < 12 ym spectral *window region are given 
in Table 2-10. 











I 
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Table 2-10. Values of the Asymmetry Factor in the Infra-red Spectrum 



AX (pm) 

8.00 

8.40 

8.60 

8.80 

9.00 

9.20 

9.40 

9.50 

10.00 

11.30 

12.10 

12.50 

1 ■ 

s’ 

.40 

.16 

.09 

.04 

.06 

.06 

.10 

.12 

.15 

.27 

.49 

.17 


It may be inferred, from the material in Table 2-10 that the truncated 
phase functions are nearly isotropic. 

Sargeant and Beckman (1973) found that the heating rate, even at the 
surface, was in error by less than .1 out of 1.9°K day ^ if they assumed 
isotropic scattering. Their error decreased strongly in relative value 
with height for a quartz aerosol over the Indian desert. 

For the sake of simplicity, we assume that the truncated phase function 
is isotropic, i.e. g* - 0. The aerosol is assumed here to be optically 
active only in the 8 - 12.5 pm region. Outside of this region, w* = 0 . 


2-7- c) Solution of the Radiative Transfer Equation 

The solution to Eq. 2-7,2 may be separated into two parts - one for the 

emitted intensity, I , and one for that scattered, I . Thus 
e s 


I - I + I 
s e 


2-7, 3a) 


M dF " l e " {1 - t ° ,)B 


2-7 ,3b) 


M dF“ r s " “ ,J o 

For later reference let us also define 


P -7-T - 1 - B 

dt o 


2-7 ,3c) 






where 


and 


J 0 ( T '.M ) -f 


p(y» m'JKt', u’My* 
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2-7, 3e) 


-1 


I Q - I e (0)' = 0) 


2-7 ,3f) 

The solutions to Eq. 2-7, 3b) and 2-7,3d) are immediate, however the equation 
for the scattered intensity must be solved by an iterative process or by 
other numerical methods. 

/ . v e have chosen to approximate the n-th order of scattering in the 
, following way 

j 

I _ n 

2-7.4 


I - 


i (T'.e) = (t\8) 

an e 


t 


where I^T'.e) is the incident emitted intensity at level r' from direction 8, 

Then the total scattered intensity is set equal to 

\ 00 

I (T’,0) - I (x',8) = I SlljSl 2-7.5 

8 n-0 e e 1-oj* 


1 


I 


l i 


The equation of transfer for the scattered intensity thus becomes 

T WW ‘ W.„>- -1 1 - I. - ^ J„(f> 2-7.6 

-1 


where J (T 
o 


■> 

-l 


M ' ) dy ' 


2-7.7 


The solutions for the upward. If, and downward, I f, scattered intensities 


are 


* 

T 


VT’.M) 


l 

I 

T ' 


“Ct— T )/p . , 

l-<J(t') J 0 (t )e W 


2- 7.8a 


♦ f -(T'-t'J/jyl ^ 

(r\- U ) - J 0 <t') e 

8 l 1-wXf) ° l w l 


2-7. 8b 
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The corresponding fluxes are 


F b (t') - 2tr 


» 

T 

"fa’ - ) - J (t')E (t '-T ')dt 1 
|,l-wXt') ° 2 


2-7. 8c 


fV> 
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- 2ir 


x 

f s£ iiH 


0 z 
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The net flux 


F n (t') - F s + (t') - f 8 I (t') 


The function J may be expressed 
o 


2J <x') - B(*' ) E (T*-T’) + 

O * 


J(t , )E 1 C|t'-T , i)dt + J(t , )E l (|T , -t , |>a 


2-7.9 


where J(t’) - (l^Ht 1 ))B(t’) 


The numerical treatment of (t * ) is accomplished in the following 


manner: 


The main contribution to the integral is from the layers very close to 

' itself, due both to the variation in J(t) and of E^(t-x). The function J(t) 

varies much more slowly with distance from x than E^(t-x). We shell therefore 

approximate the integrals in J q (x) by putting J(t) equal to J(x) and removing 

it from underneath the integral sign. Making uBe of the relations 
<*Mx) -6x 

, and 2E2(x) S 6e , where 6 • 1.67, the diffueivity 



/ 


factor, one then gets for J Q (x) , 

B 6 -€(t '-T ') -4(x*-x') 

J Q (t) = -2-e + (l-J(x’))B(x’)[l - | (e 4* 


+e )] 2-7.10 


£ 

i \ 


The first term represents the surface intensity B , singly scattered, 
reaching level x, the second represents the intensities reaching level X 
from the atmosphere below and above level X, respectively. 

2-7-d) Simple Approximation to the Scattering Term 

/ 

An even simpler approximation to the scattering source function term. 


,w’J (x'),:may be based on an analysis of the ratio, R, between the scattering 

; ° j ; 

.source function term, and the source function for emission, (1 - u* (x'))B(x ') . 

i i l 

i . ■ r . I 

( where: ■ j 

{ t ig 

-&( t' t‘ 1 B 1 _5t ' 

R(x') - W '(x'}{ 1 + “ le ^ " (“5 D-e ] } 2-7 .11 

2 E (l-w’Cx ’ ) ) 


On substituting representative values for the parameters. 


lix'i3 


one finds that 


where 


.1 < u’Cx')! .8 


.5B < B(t') < B , 


R(x') - gCx'Juftx') 


.8 < g(x') < 2.3 


and g(t') is a function of x' that may be empirically determined. 


The complete equation of transfer may be approximated by 



46 


2-7-e) The Integrations with Respect tc Height and Wave-Number 


Several schemes were tested for use in the integration with respect to 
height. It was found that an integration in the optical depth, T, in 26 steps 
of equal intervals using Simpson's rule gave numerical accuracies of the order 
of .1 percent in the fluxes. The fluxes were then interpolated parabolicallv 
in order to arrive at values for the fluxes at fixed pressures from the surface 
in 40 mb steps. 

The spectral cooling rates at wave number, v, are evaluated by use of the 
following equation 


3T v 

8.64 * 10 


4 


& 


c P 
P o 


3F(tI ?) ) 
n,v 

3(P/P o ) 


8.34 x 10 


8F n, *"<>•» . -! 

mi y “ 

O 
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The flux divergences are computed numerically. 

The integration over the spectrum between 800 and 1250 cm \ the spectral 
region where we introduce the desert aerosol, is done numerically using the 
trapezoidal rule. The results of this integration are then added to those of 
the non-gaseous atmosphere in other spectral ranges, to obtain the total 
heating rate. 

2-7-£) Comparison of Infra-red Heating Rates with Previous Comrutatlons 

We now test our approach by comparing our fluxes and heating-rates 
for a quartz aerosol to those found by Sargeant and Beckman (1973). The 
latter evaluated theoretically the i.r. heating due to a quartz aerosol 
embedded in a tropical desert atmosphere over Rajasthan, India. Temperature, 
humidity and aeroBol profiles were all measured by Peterson (1968) on April 
24, 1966 over New Delhi as was the total mass. 


The comparison is shown for the heating rates in Figure 2-3. The differ- 
ence in the heating rates for the clear atmosphere between the previous and 
present results is due almost entirely to the different models for spectral 
gaceous absorption used. In the previous study of the effect of quartz 
aerosols it was assumed that in regions of weak gaseous absorption, one could 
neglect the latter compared to the effect of the quartz aerosol. A comparison 
of the aerosol and gas optical depths as given in Table 2-11, shows that one 
should not neglect the gaseous absorption in the "window" regions. 

In the present study, we did not neglect the gaseous absorption and most 
importantly cook into account the so-called e-type absorption, supposedly due 
to dimers, see Cox(1973). This leads to an increase in cooling near the 
surface and Co the relative maximum between 950 and 850 mb. This change in 
the heating-rate profile is also reflected in the heating of our dusty atmo- 
sphere. The relative minimum in the cooling near 800 mb in the dot-dashed 
curve is due in part to the rapid decay of the effect of e-type absorption 
of the water-vapor and in part to the presence of the aerosol. The profile 
of the composite albedo for single-scattering is very different from that 
of the aerosol alone which is constant with height. These are the reasons 
that the previous authors could not reproduce this behavior of the heating 
rate profile, which is shown by the profiles measured by Bryson et. al.(1964) t 
Peterson (1968), and by de Luisi (1975). 

At the higher levels, where the effect of water-vapor abosrption is 
small, our results over-lap those of Sargeant ar.d Beckman and both merge 
with the cooling of a clear atmosphere above 400 mb. 

When we neglect the gaseous absorptions in the 8.0 ym < X < 12.5 yo 
region, we get for quartz a curve (dashed in Figure 2-3), which is very 
close to that of the previous study. This shows that under the same assump- 
tions, our method gives the same result as that c i Sargeant and Beckman. 



, HEATING RATE ( A day * 

1 l 

i 

Figure ‘‘ 2 - 3 Comparison of aerosols on thetrmal infra-red heating rates found 
by Sargeant and Beckman (1973). ( ) gaseous absorbers only. (• • •) gas- 

i 

eous absorption, with e-type water vapor absorption. ( ) quartz fros 

Sargeant and Beckman. ( ) quartz present*, analysis, no gaseous absorp- 
tion for 8.0 pa < X < 12.5 pm Q (— •— ) desert aerosol: BF • 1. 

Table 2-11. The Optical Eepth/i at the Surface of a "Normal" Atmosphere 
Containing Desert Aerosols and of the Caseous Absorbers. 


Wave length 

Aerosol 

Gases 

linn) 

8.0 

.34 

.38 

8.6 

.29 

.30 

8.8 

.32 

.28 

9.2 

.29 

.24 

9.4 

.28 

.25 

9.6 

.30 

.27 

9.8 

.28 

.29 

10.2 

.22 

.24 

10.8 

.17 

.26 

11.4 

.15 

.28 

11.6 

.12 

.29 

12.0 

.11 

.31 
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The remaining smaller differences between the two models for quartz are 
due to our different treatment of the forward peak. The previous study as- 
signed the area under the phase- function between coa0 a .75 and cos8 ■ 1 to 

2 

the forward peak, while we used g for this area. The areas at a wavelength 
of 10 pm are .40 In the present case and .46 previously. Secondly, we 
assign a best-fit scale height “ 2.8 kilometers and a cut-off at 400 mb 
to the height-distribution instead of the actual distribution. Thirdly, we 
assume that the aerosol effect may be neglected outside of the spectral 
window of water-vapor between 8.0 pm :£ X < 12.5 pm. 

In summary then it may be said that our method of computation compares 
well with the previously described method when applied to the same data. 

2-7-g) Infra-red Heating Rates < n a Desert Aerosol Layer Aloft Over the Ocean 

We have modelled the effect on thermal infra-red heating rates of a 
desert aerosol overlying a tropical ocean in a layer between 2-5 km. This 
Is close to the situation existing from late spring to late summer over the 
tropical Atlantic Ocean. The amounts or masses cf aerosol that wc Iiav?. used 
are close to those which have been observed for the light and noiral cases 
(BF - .1 - 30 Pg/m 3 8 fl (.5 pm) - .019; BF - I. • 299 Pg/m 3 , 3 a <.5 pm) - .193) 
in the Western Atlantic. The heavy desert dust (BF «* 10.) may occasionally 
be found over the eastern tropical Atlantic off the coast of Africa. The 
heating rates caused by this latter aerosol constitute an upper limit on 
the possible effects of the desert aerosol. 

Figure 2-4 shows the location of the layer between 2 and 5 km. The 
full line and the dotted line show the cooling in a clear tropical atmosphere 
with and without e-type absorption, showing the possible importance of 
this effect in the lowest 200 mb of this humid atmosphere. A chin aerosol 




(BF - .1) will reduce cooling by about .1 to .2*K day 1 below the aerosol 
layer and in its lower part, and by about the same amount of additional cooling 
in the upper-part of the layer. The normal case (BF *1.) will reduce cooling 
up to 1°K day - below 750 mb and increase it up to .75°K day in the upper 
levels of the dust layer. A heavy dust fall (BF «* 10) will have a very large 
effect on the heating rate, leading to strong heating around 800 nib and large 
cooling in the upper part of the dust layer. 

In every case, the resulting cooling profile would lead to a 
stabilization. 


1 


4 4 -« -l -1 -I O’ I > > * » • 1 * 

mcatim* natc 

Figure 2-4 Infra-red heating rates for an aerosol cloud over the ocean 
Desert aerosol aloft between 2 and 5 km indicated by shading. Tropical 
atmosphere, McClatchey (1S71). ( ) gaseous absorbers only. (•*•) gas- 

eous absorbers only without e-type absorption for water vapor. ( — •••• — •) 

BF » 0.1 light desert aerosol. ( ). BF - 1.0 normal desert aerosol. 

( — ) BF - 10 heavy desert aerosol. 
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2-8 The Total, Solar and Infra-red, Heating Rates 


Combining the considerations presented in sections 2-6 and 2-7 yields 
the total heating rate due to solar and infra-red radiation. These computa- 
tions are applied to vertical aerosol distributions described by the parameters 
cited in Table 2-6. 

2-8-a) Desert Aerosol Over Land 

Figure 2-5 shows the total daily heating, evaluated for three 
amounts (BF ■ .1, 1., 10.) of desert aerosol over New Delhi using the 

, i \ 

atmospheric land aerosol profiles for April 24, 1966 measured in the evening 

i ,j ' i 

by Peterson, |(1968) , and for a surface albedo of .25. Most of the effect of 

i ! : 1 

the aerosol is to be found in the lowest 3 kma of the atmosphere in the 

4 , I 

case cf b light and a normal desert aerosol. The main result of the presence 
of a light aerosol is an increase in cooling between 2 and 3 km of the order 
of ,5°K day *. A normal aerosol will strongly reduce the cooling between 
.5 and 3 kms by 1.5®K day * and slightly increase it by .25°K day near the 
surface. A very large amount of dust will lead to strongly Increased cooling 
near the surface, strong heating between 1 and 4 kms and reduced cooling all 

the way to 10 ’ms, leading to a stabilization process. 

The average cooling per day for the atmospheric column between 0-10 km 

is -1.68°K day * in the clear atmosphere and -1.86, -1.41 and -.18®K day 
for the case of a light, a normal and a heavy aerosol, respectively. The 
results show that it is danger*, us to make generalized predictions about the 
effect of aerosols on atmospheric heating. One must take into account not 
only the surface albedo and the amount of aerosol, but also its vertical 
distribution and the amounts of the absorbing gases. 
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HEATING RATE Ckdoy' 1 ) 
Figure 2-5 Total dally heating rate profile over a desert. 

A « 0.25, ■ 2.8 km ( •) clear atmosphere. (••••) BF ■ 0 

BF - 1.0, normal aerosol. ( ) BF - 10, heavy aerosol. 


New Delhi, April 24, 
.1, light aerosol ( 


1966 

) 


2-8-b ) Desert Aerosol Over Ocean 

Figure 2-6 shows the results of similar '’alculations for the 
case of a desert aerosol aloft between 2 and 5 kins overlying a tropical ocean 

with a surface albedo of .05. Again, the effect of a light (BF ■ .1) aerosol 
is small (.15 - 2°K day *) but systematic. A normal desert aerosol has a 
much larger effect. The daily cooling is reduced by .5°K day ^ in the lowest 
km, there is an actual heating of .9°K day ^ at the bottom of the aerosol 
layer and a broad increase of cooling of about .4°K day ^ in the upper part 
of the aerosol layer. 
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The daily cooling from the ocean surface to the bottom of the aerosol 
layer ia -1.45°K day * in the clear case, -1.36*K day * for a light aerosol 
and -.78*1^ day * for the normal case* Inside the dust layer itself, the 
change in cooling is small -.91, -.83, -.89°K day * for the light, normal 
and heavy aerosol loads respectively. The important effect thus seems to 
be the reduction in cooling below the dust layer and in the latter* a lower 
part. 



DAILY HEATING. Pk doy'*) 

Figure 2-6 Total daily heating rate profile over the 
ocean. Ar. aerosol layer between 2 and 5 km is designated 

by shading (••••) clear atmosphere, ( ) Br - 0.1 light 

aerosol. ( ) BF = 1.0 normal aerosol. 

It is of great interest to measure the aerosol concentration in those 
regions of the atmosphere where one might expect heating rates that are 
different by . 5°K day ^ from those of the clear atmos'pnere. Such measure- 
ments can in principle be achieved by a Shuttle-borne lidar. 
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3-0 Sensitivity Requirements for the Determination of Desert Aerosols from a 
Shuttle-borne l<da r 

We will address ourselves to the determination of the accuracy with which 
the aerosol extinction coefficient must be measured and the vertical resolution 
with which such measurements must be determined. It may be recalled from the 
discussion in section 2-1, b) that G.C.M.s are Insensitive to heating rate 
perturbations less than 0.5°K day The aerosol concentration must there- 
fore be determined so that the computed heating rate profiles are defined 
within those tolerances. The Shuttle-borne lidar will measure profiles of the 
aerosol extinction coefficient; we assume that such measurement* are obtained 
at A « 0.5 pm and that these profiles are uniquely related to aerosol con- 
centration as predicted by the aerosol model defined in Ch. 1. 

The aerosol extinction at A ** 0.5 pm and the accuracy required to stay 
within the heating rate tolerances were computed as a function of altitude. 
Although 0.5°K day * is a tolerable heating rate uncertainty for a G.C.M., 
we also present results for 0.2°K day 1 perturbation. We will show in Ch. 4 
that the p ‘sent state of technology can yield extinction coefficient profiles 
which are consistent with heating races defined within 0.5°K day The 
advance of technology may render it feasible to define heating rates within 
0.2°K day * by the time that the Shuttle is launched; the corresponding 
extinction coefficients are therefore also presented for that case. 

Calculated heating rate profiles are presented as a function of altitude 
for various aerosol models characterized by parameters cited in Table 2-5. 
Results are presented for aerosol profiles characterized by scaling factors 
BF ■ 0, 0.1, 1, 5 and 10. It la noteworthy that BF ■ 1 and 5 are similar to 
the clear and hazy tropical model atmospheres of McClatchey et al. (1971). 
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These models yield visible ranges at the surface of 22 and A. 5 km. respectively. 
Although results for total and solar heating rates are presented for various 
aerosol models in Tables 3-1 and 3-2, the analysis of Ch. A will be confined 
to the BF «* 1 and 5 cases. 

We carried out another test in order to determine whether a resolution 
of 1 km in the determination of the extinction coefficient will change the 
atmospheric heating rate, by more than . 2 e K day. An example of our results 

I i 

’ t < 

is shown in Fig. 3-1. Part (a) of the Figure shows the heating rates in the 

' ■ 'I : ■ { 

solar visible spectrum (.3 pm < X < .8 un) for the caBee of a light, a normal 

l ■ 'i i 

l ■ t 

land a heavy desert aerosol (BF ■ .1, 1 and 10 respectively). Part (b) of the 

same Figure shows the error in degrees per day when the average heating rate 

\ 

in successive layers of 1 km thickness is applied at the center of each layer. 
The error, as typified by the difference between the average heating rate in 

\ 

a one km layer and the actual heating rate at the mid-point of the layer, is 
larger than 2 percent only very near to the surface for a normal and a heavy 
desert aerosol. The largest absolute value of the error is only .055°K day \ 

It seems therefore that a resolution of one km in the lidar measurements is 
adequate for the determination of the extinction coefficient with an accuracy 
that is required by the G.C.M.s. 
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Table Ch«nge in R^ at 1 • 0.5 t«a to Produce a Change of .2*t day 

In Total Dolly Meriting In a busty Ateoaphcra 



BP - 0 


BP - .1 

BF - 1 


tf - 5 

BP - 10 

Height 

|ko) 

>T 

at 

8 . 

|Af>J 

dT 

at 

8 

a 

|M«I 

3T 

at 

6 . 

l M .l 


3T 

at 


lAB.I 

AT 

at 

*. 

IM.I 

0 

- 

.87 

0 

.19 


- .83 

.020 

.05 

' 1.09 

.19 

.08 


-4.55 

.97 

.13 

-6.03 

1.91 

1.93 

1 

- 

1.4d 

0 

.20 


-1.31 

.014 

.14 

- .57 

.14 

.06 


- .37 

.70 

.43 

- .21 

2.33 

. 89 

2 

- 

.98 

0 

.01 


-1.25 

.009 

.01 

- .03 

.09 

.02 


♦ 1.73 

.45 

.09 

♦2.53 

.94 

.19 

i 

- 

.46 

0 

.10 


- .48 

.007 

.06 

- .63 

.07 

.06 


« .35 

.35 

.13 

♦ .74 

.66 

.17 

4 

- 

.53 

0 

.05 


- .5* 

-C03 

.05 

- .70 

.05 

.04 


- .30 

.15 

.05 

♦ .32 

.46 

.14 

5 

- 

.66 

0 

.07 


- .62 

-C03 

.03 

- .63 

.03 

.03 


- .33 

.10 

.11 

- .12 

.33 

.20 

6 

- 

.90 

0 

.04 


- .89 

.002 

.02 

- .82 

.02 

.03 




- .68 

.08 

.05 

- .36 

.23 

.14 


Itotti 

Ac«oa: Hew Delhi April 24, 1966 

Surface Albedo - .25 

Scale Height of Desert Aeroaol • 2.8 ka 
Top of Aarosol Layer at 400 ob or 6.7 ka 


Iable 3-1-b) Change in it X ■ 0.5 la to Produce a Change of .5*1 day 
In Totel Dally Heating In a Dusty Atmosphere. 



Ataoe: New D.-lhl Ap**l 2* , 1966 

Surface olheJo • .25 

Scale height of desert aeroaol *s 2. 8 ka 
Top of aerosol layer at 400 th ot U la 
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Figure 3-1: (a) Heating rates due to desert aerosol 

in the solar spectral region .3 - .8 yra. T* <= .031 - 
BF - .1; T* - .305 - BF » 1.0; T* - 3.052 - BF - 10. 
Surface albedo * .25; Cosine of solar zenith angle ■ . 
Aerosol Scale Height » 2.8 km. 

(b) Difference between heating-rates 
averaged over layers of 1 km thickness and the heating 
rates at the mid-point of the layer in degrees per day 
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4-0 Lldar Methods to Measure Aerosol Profiles - Introduction 


The criteria required to define the spatial distributions of tropo- 
spheric aerosols for input into atmospheric general circulation models were 
identified in Tables 3-1 and 3-2. We will now describe several schemes in 
which lldar configurations are used to measure the spatial distribution of 
aerosol optical depth. Error analyses and model atmospheric calculations 
are presented in an effort to ascertain whether these schemes can provide 
spatial distributions of aerosol optical depth from Shuttle missions. 


4-1 General Discussion of Expected Lldar Signals 
4-1-a) Generalized Lldar Equation 

The lldar equation for the number of photons collected by the receiver 
telescope due to single backscattering from a horizontally homogeneous atmo- 
sphere in the presence of J distinct extinction components 0^ (J - 1, J) at 
an angle 0 to the vertical is 


2 ( x ) 

c dc 


Vr V U 1 


hcO 


4tt 


sec z 0*x 2 


expj -2 sec0 


J \ 4 - 1,1 

) 


T 

J 1 




J* 


where x - equivalent vertical path length as measured from the lldar platform 
c - speed of light 
E & - energy of optical output 

hcO - photon energy, where 0 * 1/A is the optical wavenumber and h is 
Planck's constant 


A “ receiver area 
r 


+ aj) * albedo for single scattering, where s^ and a^ denote 


extinction due to scattering and absorption, respectively 


IP, (it) 


4n 


normalized backscattering phase function and 


" S J + a j ” C ^ e extincCion coefficient of the jth component. 


IPjOO and {3j are general functions of equivalent vertical path length x. 


The geometry of £q. 4-1,1 is shown in Figure 4-1. 



HORIZONTAL RANGE 

Figure 4-1 Lidar equation geometry for the case 
of a horizontally homogeneous atmosphere. 


A lidar detection system integrates this signal over a finite time 
A/ 

interval At, so that the number of singly scattered photons collected by the 


receiver telescope from the equivalent vertical range x to x + u is 
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x o + u 


*j<V u) 


E A 
o r 


, dN (x) 

— — -J sec0 dx 

c dt 


hco secQ 


+ u 

a ip. (it) e.w. 


4it 


( 

*0 

J 

il -28ec0 

z e.,dx’i 
J'-l J 


4-1,2 


dx 


cAt 


where u “ 2 ~ sec 6 * ^* e nusa ^ er photons collected by the receiver telescope 

dN (x) 

— — is approximately constant in the range interval u if u is sufficiently 


small. Equation 4-1, 2 then reduces to the more common form 

I 

E A u IP. (it) B.w. 


o i 


V X o* hcO sec8 4 tt 




exp \ -2sec0 


4-1, 3 


Z 0, ,dx* 


This chapter presents investigations of lidar returns from an atmosphere 
composed of Rayleigh (molecular) and Mie (aerosol) scatterers, which are 
referred to by the subscripts m and a, respectively. 


4-1-b) Lidar Returns from a Model Atmosphere 

Lidar returns are presented for a downward viewing lidar using the Clear 

and Hazy Tropical atmospheres of McClatchey et al. (1971). Example height 

profiles of B and 8 for two Tropical model atmospheres at 4880& are 
m a 

plotted in Figure 4-2. The height profiles of for the cases of Clear 
and Hazy model atmospheres differ only in the lowest 5 km. 

Figures 4-3a and 4-3b present lidar returns from the McClatchey Clear 
and Hazy Tropical model atmospheres, respectively, for a lidar platform height 
h Q of 185 km. The lidar platform is located at a height of 185 km to simulate 
a Shuttle borne lidar system. The returns are given in units of photons per 
milll joule output at the wavelengths 3371, 4880 and 6238&. The returns are 


McCLATCHEY TROPICAL 
ATMOSPHERE X=48802 
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Figure 4-3a Lidar returns N and N from the Clear Tropical 
j o a r 

I model atmosphere of McClatchey et al. (1971) at 3371, 

ij 4880, and 6238 X. The lidar platform ie located at a 

j I . 1 

l height of 185 km to simulate a Shuttle-borne lidar 
! ' ■' i 2 

system. The receiver area ■ 1.0 m . 



Figure 4-3b Lidar rer srne and N a from the Hazy Tropical 

model atmosphere of McClatchey et al. (1971) at 3371, 

4880, and 6238 X. The lidar platform is located at a 

height of 185 km to simulate a Shuttle-borne lidar 

2 

system. The receiver area A " 1.0 a . 
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separated into the quantities N and N , which represent the number of photons 
collected by the receiver telescope due to single backscattering from molecules 
and aerosols, respectively. The returns of Figures 4-3a and 4-3b were calculated 
using a vertical range interval u * 500m and a receiver area A r " 1.0 m . Fig- 
ures 4-4a and 4-4b present lidar returns in the same format but with a platform 

height h of 10 km to characterize returns measured from an aircraft. The 
o 

returns of Figures 4-4 were calculated using a vertical range interval u ■ 200m 

2 

and a receiver area A - 0.1 m . The optical properties IP (ft) and to used 

r a a 

for these calculations are given in Table 4-1. 



Table 4-1 

The Model Atmosphere Properties 
IP (ft) and 0 ) tor Figures 4-3 and 4-4 

A A 


Wavelength 

ip (ft) 

a 

0) 

A 

3371 X 

0.40 

1.0 

4880 X 

0.40 

1.0 

6238 X 

0.35 

.98 


j 

4-2 Signal and Noise Considerations j 

■! 

4-2-a) Safety Standards 1 

1 

The Department of the Army and Navy have set guidelines for the ; 

Maximum Permissible Exposure (MPE) to laser radiation, see Pressley (1971). 

The Army and Navy standard maximal exposure level which is not expected to 
cause detectable ocular injury (measured at the cornea) is about 1 x 10 

I 

j cm ^ for a Q-switched laser at all wavelengths between 0.4 and 1.5 lim. 

A safety factor of 2 is recommended for use in field evaluation. 



67 



v 

h 


5 

i 


!' 

r 



Figure 4-4a Lldar returns N and N from the Clear Tropical 

model atmosphere of McClatchey et al. (1971) at 3371, 4880, 

and 6238 X. The lidar platform is located at a height of 

10 km to simulate an aircraft-borne lidar system. The 

2 

receiver area A «■ 0.1 m . 



Figure 4-4b Lidar returns N and N from the Hazy Tropical 

m a 

model atmosphere of McClatchey et al. (1971) at 3371, 

4880, and 6238 X.' The lidar platform is located at a height 

of 10 km to simulate an aircraft-borne lidar eyatea. The 

2 

receiver area A ■ 0.1 m . 




^As stated in Moses (1970), the diameter of the normal human adult 
pupil under ordirary room light levels is between 3 and 4 mm. It is further 
noted that retinal sensitivity under daytime exposure levels is consider- 
ably less than for a normal dark adapted eye. We therefore infer that the 
ANSI MPE values for visible and near infrared radiation are conservative 
for daytime laser exposure in the outdoor environment. However, the ANSI 
Standard section 8.5.1 states that no corrections are to be made for the 
intrabeam viewing MPE when the limiting aperture (pupil) size differs 
from . 7 mm. 
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4-2-b) The Impact of Safety Con ;i 1 derations on the Lldar Field of View 

The following calculations provide estimates of the minimum allowable 

laser transmitter beam divergence for 0.4 Mm < A < 1.4 pm. The ANSI (1973) 

-7 —2 

maximum permissible exposure per pulse of 5 x 10 joule cm at the Earth's 
surface Is chosen as the eye safe value for lidar operation in an outdoor 
environment . 

The instantaneous radiant exposure at a distance from a laser trans- 
mitting energy through the atmosphere will be determined by 

1) atmospheric attenuation due to scattering and absorption, 

2) the distribution of optical energy within the transmitted beam, , 
and 3) atmospheric scintillation, whi-.h cause random and possibly hazard- 
ous radiant exposure maxima. 

Atmospheric attenuation can be neglected for the purpose of obtaining a worst 

case estimate of surface radiant exposures. 

The distribution of optical energy within the transmitted beam is not 

uniform over the entire field of view ft Q . If the emitted energy were 

distributed uniformly over the beam area, the mean radiant exposure per 

pulse J at range x would be 
max o 


E JrV 2 
o o o 


4-2,1 


where E is the total energy emitted per pulse within the solid angle ft . 
° o 

The Gaussian is a more realistic model of the energy distribution within 
the beam area, see Figure 4-5, The Gaussian energy distribution E(C), 
where C is the half-angle from the transmitter optic axis, is specified 
by the requirement that the fuil beam divergence is that solid angle 




A- „J 



Figure 4-5 A comparison of the uniform and Gaussian j 

models for the angular distribution of output laser < <" 

energy. The shaded area represents that portion of \ 

the Gaussian distribution (86.5%) which falls within ’ 

the full angle beam divergence. ! / 


which contains ^86.5% of the transmit' ed energy. E(C) 1b then given by 
E a 

E(C) - "J- exp(-aC 2 ) 4-2,2 

, 2tt 1 , 

for a ■ jj- • 4-2,3 

o 

The energy per unit solid angle at the center of this distribution is 



i 


i 


i 

j 



A comparison of Eqs. 4-2,1 and 4-2,4 shows that the more realistic 
Gaussian energy distribution gives a peak mean radiant exposure which is 
twice as great as that due to the uniform energy distribution. 
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The Instantaneous radiant exposure will depart significantly from its 
mean value due to scintillation of the beam during transmission through the 
atmosphere. A useful estimate of the probabiliti.es for encountering 
hazardous radiant exposures due to random beam fluctuations is given by 
Dabbert and Johnson (197i). Their experimental work determined the effect 
of range, laser wavelength, atmospheric thermal turbulence intensity and 

I path geometry on the magnitude of laser scintillation. Figure 38 of 

i Dabbert and Johnson, ■ which is itself a conservative estimate, shows that 

' ' I -2 

! the instantaneous probability for eyr. damage will be less than 10 when 

. ij \ 

I the mean radiant exposure is an order of magnitude less than the eye safe 

i ; ' i 

value. The results of Dabbert and Johnson are given in Table 4-3. We con- 
clude that a radiant .exposure which is a factor of ten below the ANSI (1973) 
MPE limit is adequate for eye safety. 



I 


l 


Table 4-3 Worst Case Instantaneous Probability for Eye Damage due to 

Atmospheric Scintillation as 

Calculated by Dabbert and Johnson (1971) . 

Ratio of eye safe to mean * 
value or radiant exposure 

Instantaneous probability for 
exceeding the eye safe value 

> 10 

< 10 -2 

25 

10‘ 3 

50 

10~ 4 

100 

ic -5 

200 

10- fc 

t'fhese results are based on a 

5 mm limiting aperture 
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We conclude from the above that the Instantaneous radiant exposure 
which is eye safe at the Earth's surface is 


J , • 2E fi * 2 S J , x 10 -1 

sfc o o o safe 


4-2,5 


Using the ANSI MPE value J flafe “ 5 x 10 7 Joule ca" 2 (See Table 4-2), we find 
that the laser transmitter beam divergence must exceed the following minima: 


fl 


o,min o 


E ^ > 20J_* c _h_ 2 = 1.2 x 10 7 star ad joule ^ for h Q ■ 185 km 4-2, 6a 


safe o 


-5 -1 

- 4 x 10 sterad joule 


for « 10 km 4-2, 6b 


The practical limit for laser transmitter beam divergence is 
* -9 

0 . * 1 x 10 a ter, which occurs at 

o,min 


E = 8.3mj for h ■ 185 km 
o o 

= 0.025mj for h ■ 10 km 
o 


4-2, 6c 
4-2, 6d 

Note that Eqs. 4-2,6 imply a safety factor which is two times safer than the 
Army and Navy field evaluation criterion given in section 4-2-a). 

Assuming that the laser output is axially symmetric, the full angle divergence 
equivalents to Eqs. 4-2,6 are 


, = .38 

full, min 


•7.1 ] 


E,, 

mrad 

for h - 185 
o 

Iqjj 

4-2, 7a 

1 joule 

Eo 

mrad 

for h “ 10 

o 

km 

4-2, 7b 

1 joule 


The practical limit for the full angle divergence of the laser output is 

= 0.036mrad 4-2, 7c 

which is independent of transmitter altitude. 


* 

^full.min 


The area per unit energy of the illuminaced "spot" on the Earth's surface is 
-1 _ 2-1 -l t - -i 


A , E„ 
surface o 


« h 2 e _1 

o,min o o 


20J safe = 4 x 1()3|: joule 


4-2, 8a 

which is independent of transmitter altitude. The practical limit for the 


4-2, 8b 
4-2, 8c 


crea of the "spot" is 




Cin h o 2 - 34 ° 2 

for h * 
o 

185 

kin 

•* 0.1 m 2 

for li - 
o 

10 

km 



The worst case daytime solar signal can be estimated by scattering all 
of the available solar radiation from a Lambertian surface which completely 


fills the receiver field of view The worst case solar signal estimate, 

given in number of solar photons collected by the receiver teleccope, is 




I X 1 AX 

N _ — (— - — cos0 )A At 
s ir hco art 


4-2,9 


where I. is the solar spectral irradiance at 1AU, AX is the receiver band 

; A , 

; I 1 <\i 

| pass at wavenumber o, 0 is the solar zenith angle. At is the range gated 


(signal averaging time, A is the albedo of the Lambertian scattering surface, 

I ‘ i 


' c is the speed of light, and h is Planck’s constant. Atmospheric attenuation 

' of solar radiation haG been neglected. A tabulation of the solar spectral 

irradiance at 1AU is given in Table 4-4. 

The backscattered signal N is estimated by subtracting an independent 

measurement of from the total lidar signal. The signal to noise ratio 

(SNR), considering statistical variations in the backscattered and solar 

signals but neglecting noise contributions due to system electronics and 

— 

optical processing, is given by N*(N + N^) . Since N <* E q and 

N a Q . * E for n > Q* , when the lidar is operated at the worst 

6 o,min o o.min - o,min 

case eye safe limit, the signal to noise ratio is effectively proportional 
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Worst case estimates for N at representative wavelengths are given in 

Table 4-5. Since ft . is proportional to E , these N estimates are given 
o,min o 8 

in units of cumbers of photons per milli joule output. The High Spectral 

Resolution Lidar system has a bandpass SX^ or 60 ^ - fiX^ * X . Values for 

<5Xp at an atmospheric temperature T ■ 250 e K, and the solar signal which ie 

passed by such a filter are also given in Table 4-5. For E < E . ft 
r J o o * o,mln 

* 

cannot be reduced beyond the practical operating limit ft . and the solar 

c ,min 

signal assumes its minimum value N . . Shuttle case values for the worst 

s, min 

fg 

case solar signal at minimum practical beam divergence ft . are also given 

o,min 

in Table 4-5. 

4-2-d) Reduction of the Solar Signal through Lidar Operation in a 
Fraunhofer Line . 

The Table 4-5 values for the solar signal can be reduced through lidar 
operation within a Fraunhofer line. A partial list of suitable Fraunhofer 
lines is given in Table 4-6. Included In this table are: 
l; peak line depths in percent of the solar continuum, as measured at the 
Earth's surface for direct solar radiation, and 
2 ) the line width measured at twice the peak line depth. 

Figure 4-6 presents the location of several of these lines with respect to 
the output of a commercially available dye laser. 
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co# 6^ * 1 


I 

0 • 1.2 x 10 Joule output \Eq. 4-2,6#) 

‘ -v 

At - 0.31 * 10“ 5 ##c - 500n) 


lb# values for * AX * and * AX * r « based on single saaples for eye m/« operation at 1 *J 

l.a.r output. For total solar .l,nal bated on an average of H eamplaa with • lour output of E (JouluJ, 

o 

multiply these values by Vt Q / 1 The following values are used for the aircraft cm: 

h - 10 hm 0 • 4 x 10~* star joule 1 witput (Eq. 4-2, 6b) 

° 2 r 
K % - 0.1 ■ 

far the aircraft cm, divide the solar signal values by 3. Lldar operation in a rreuniwfer line #ay 
raduca the solar signal by an order of magnitude. 

(2) From Table V ©f HA5A SP-6005. Kay, 1971. £ I^dX - 1351 watts u~ 7 . 

(1) The Doppler width (FVKM) (or photons which ere backscattered by molecule# la 


1351 watts u ~ 7 . 


tk - * 

I) c Y _ 


where k Is dolttoenn's constant. ■ Is the areo molecular *■•« end T Is the atmoepherlc temperature 
(1 was chosen lo be 250*1). Note thee Ac^ *• 61^ • X 7 . 

s 

(4) Toe its* fitter bean divergence sod receiver ileid of vi«w have the practical Halt ^ • 0 f 

• 5 x Iff* star. This limit is reschod when f i 42 sj for the Shuttle cess, and whan E < .13 sj fot the 

o o 

aircraft case. For the aircraft case, divide the R # tt4n values by 10. 


Table 4-5 Typical Veluea for the Iwicr Si^iaJ at Representative Wavelength* - Shuttle Case* 




Coucerin 2 
USU-4750 X 


Rhod amine 6C 
4500-6500 X 


(oS» 

!«~*l 

0 <*> 
«,*io 

(photons) 

iw.2 

.125 

i) 

16.9 

.093 

46 

20.3 

.086 

50 

21.6 

.082 

55 

23.1 

.076 

58 

26.2 

.067 

71 

29.2 

1 

.061 

76 

i given in unite 


(1) Th. aol.r «l,o.l v.lu*. «c« calculated for aloloua bu. divergent. , ng or. (Irm 1. unit* 

I j 

of photon. P«r olllljoul. .'of tranaalcted laa.r eoergr- to. tolar .i,o.l 1. that mjabor of aol»r photon* 
coll.c t ad bp th. r.cl.ar t.l.aeop., Th. calculation* uaa tha fellowin, Shuttle uu ualua.i 
















78 


A number of observers have noted that the Fraunhofer lines seen in 
atmospheric diffuse radiation are relatively less deep than in direct solar 
radiation. The filling-in of the Fraunhofer lines is not fully understood t 
and the effect is given the generic label "Ring effect" for convenience, see 
Grainger and Ring (1962). The Ring effect is thought to be primarily due to 
flourescence of air or airborn aerosols (or both). A portion of the effect 
could also be explained by rotational Raman scattering by Nj and 0^. 

As cited by Barmore (1975), the magnitude of the effect is 
reported to be a few percent of the adjacent 3olar continuum: 

Grainger and Ring (1962) find the H Fraunhofer line 
at 3968 X filled in to the extent that the Ring effect 
would be approximately 5% of the adjacent continuum. 

Noxon and Goody (1965) note that for observations taken 
near noon, the Ring effect is about 3% at 4300 X and 
1.5Z at 6563 Xf 

The Ring effect shows variations with wavenumber, solar zenith angle, the day 
and time of observation, and the surface albedo, see Hunten (1970). To our 
knowledge, measurements of the Ring effect with downward viewing spectrometers 
have not been conducted. In the absence of measurements, we assume that the 
Ring effect is comparable in magnitude in both the upward and downward direc- 
tions. We conclude that the solar background radiation can be reduced by 
at most a factor of 10 through the operation of lidars within . Fraunhofer lines. 
4-3 High Spectral Resolution Lldar 
4-3-a) Introduction 

The conventional elastic backscatter lidar measures 

N (x ,u) + N (x »u). 
m o a o 

^Barmore (1975), pp. 1491-1492. 
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A high spectral resolution lidar system, or HSRL, which is capable of 
measuring N (x ,u) and H (x ,u) separately is investigated. 

The spectral distribution of 
backscattered photons near the out- 
put laser wavenumber is due to 
1) spectral distribution of trans- § 
mitted laser light, 2) Doppler 
broadening of the Rayleigh scattered 

' I 

photons due to thermal molecular 
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motion, 3) Brillouin scattering from 

. II I 

density inhomogeneities (sound waves) 

I , - s 

and 4) rotational Raman scattering 
from polyatomic molecules. j 


Figure 4-7 

Spectrum of echoes from air. 
Contributions of molecules from 
aerosols can be separated, from 
liocco et al. (1971). 


Figure 4-7 shows the elastically scattered aerosol return maximum, and 
the Doppler broadened return from atmospheric molecules. This data was obtained 
by Fiocco et al. (1971), who measured the spectral distribution of laser echoes 
about the output laser frequency using an Ar + laser (X » 4880 8). The Ar + 
laser output was stable to i 50 MHz (± .0016 cm "*■) at 0.5 watts. The piezo- 
electrically scanned Fabry-Perot interferometer had a free spectral range Q of 
4.8 C.Kz (0.16 era *) and a finesse n f ot about 12,^ see Mack et al. (1963). 


The resolving width (full width at half maximum, or FWHM) is therefore 
6o “ Q/n^ = *013 cm “, i.e. about 8 times the output wavenumber stability. 


4-3-b) Characterization of atmospheric Scattering of Laser Light 


The distribution in one dimension of molecular velocity for a gas at 
atmospheric pressures is approximately Maxwellian, such that 


dn 
m 

dv 


n 

"m V 2nkT 


exp 


r 521 
1 ~ 2kT 


where n ■ molecular number density 

m 

m “ mean molecular mass 
• v « molecular velocity 
k “ the Boltzmann constant 
and T ■ atmospheric temeprature. 


4-3,1 


According to Fiocco and DeWolf (19b8), Eq. 4-3,1 holds throughout the atmosphere 
as an approximation with about a 4% correction needed at surface molecular 
densities due to Brillouin scattering. The Doppler shift for backscattering is 


x 

o 

o 


(e - o o ) 



c 


such that 


dv [j 1 . c 

do " 2 o 

o 


4-3,2 


4-3,3 


where c is the speed of light, fT is the backscattered signal wavenumber and a 

o 

is the wavenumber of maximum emission by the laser. The spectral distribution 
of photons scattered by molecules is therefore 


me 


1 dN m 1 

\ Jf 

o 


exp 


- (0 - 0 ) 2 

8o 2 kT 0 


4-3,4 


where is the total number of photons backscattered by molecules. 
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For purposes of this discussion, the laser output is assumed to 
follow a gauss Ian distribution in wavenumber, such that the spectral distri- 
bution of photons elastically scattered by aerosols is 


N 


dtt 
a 

do 


fa 


exp 


- b (a - a o )‘ 


i-3,5 


where is the total number of photons backscattered by aerosols, and the 
full spectral width of the laser is Sj Q - . The definition Eq. 4-3,5 

assumes that spectral dispersion due to scattering by aerosols is character- 
ized by a Dirac delta function. Note that by definition 
.00 


dN 
m 

do 


do 


dN 

a 

do 


do 


4-3,6 


' l dH m » dN -* 

whcre N~ ~dO 3nd N 

m a 


are defined by Eqs. 4-3,4 and 4-3,5, respectively 


Maximal tropospheric aerosol vertical velocities may typically be 
*vin sec This corresponds to a Doppler shift which is approximately 1/100 
of the laser output resolution of a typical HSRL* see Table 4-9. Therefore, 
the Dirac delta function approximation for the spectral dispersion due to 
scattering by aerosols is good for the case of a nadir pointing HSRL. Caution 
must be exercised, however, that the HSRL does not operate at some non-zero 
angle to the nadir, or that there is a sizable vertical velocity in the plat- 
form trajectory. 

4-3-c) Brlllouin Scattering 

The spectral intensity of light scattered by a gas assumes the charac- 
teristic Doppler gausaian shape whenever the interactions between gas molecules 
can be disregarded. In this case (pressure « 1 atm), the molecular velocity 
distribution is dominated by the diffusion of thermal energy, giving rise to 
a central component in the scattered energy spectrum . As the gas pressure is 
increased, molecular velocity correlations arise from molecular interactions 
which spectrally redistribute the scattered energy. In this case (pressure > 1 


atm) , 


the propagation of pressure waves through the scattering medium gives rise to a 
symmetrically shifted Brillouin doublet. 

Yip and Nelkin (1964) have calculated the spectral response of a fluid 
model in which particle number, momentum and energy is conserved. The shape 
of the fluid model response is e function of the parameter 


a / 2kT 

2irK , „ 
m 


-*5 


4-3,7 


; which represents the ratio of the fluctuation wavelength to the collision mean 

free path. K is the wavenumber associated with the momentum transfer in the 

! 3 

scattering j 


K*2o o sin (6/2) 


4-3, 8a 


where O q is the wavenumber of the incident radiation and 8 is the scattering 
angle. The effective collision frequency a can be approximated by 


5 „ 2k 2 T 


-3 


4 -3,8b 


where p is molecular density (cm ) and is thermal conductivity of the gas 
(erg cm ^ sec * °K *) . 

In terms of the parameter y, the pressure fluctuations are collision 
dominated hydrodynamic processes when y >> 1. For y « 1. collisions are no 
longer important and the fluid behaves as a system cf independent particles. 

The results of Yip and Nelkin for the spectral intensity of light scattering 
from an ideal gas are presented in Figure 4-8. Values for the parameter y 
at several altitudes for the McClatchey et al. (1971) Tropical model atmosphere 


are given in Table 4-7. 



Figure 4-8. 

The spectral intensity of light 
scattered from an ideal gae for 
various values of the parameter 
y. The abscissa is expressed 
in terms of a normalized wave- 
number shift 6 o/6ct d , where 
600 is the wavenumber FWHM for 
pure Doppler scattering. The 
parameter y 0.3 at, surface 
atmospheric densities. 



Table 4-7 Values for the Parameter y at Several Altitudes for the 
Model Tropical Atmosphere of McClatchey et al. (1971), X • 4880X 

altitude 

[km] 

pressure 

[mb] 

temperature 

[•K] 

► ( 2 > 
parameter y 

0 

1013 

300 

.273 

2 

805 

288 

.229 

5 

559 

270 

.173 

10 

28S 

237 

.105 

20 

57 

207 

.025 


^To compute the parameter y at other wavelengths, note that Kqs. 4-3,7 

and 4-3, 8 a give y X. 
io\ 

The thermal conductivity for air can be approximated by k ^ • mpCpD, where 

c is the specific heat at constant pressure and the self-dif fupion 
P 

2 -1 

coefficient D = 0.18 cm sec for air at 0®C and 1000 mb . D varies with 

1.81 —1 

temperature and pressure by the relation D “ T p , 

Meteorological Tables (1949), Table 113. 


see Smithsonian 
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Restating Fiocco and DeWolf *1968), Yip and Welkin's results indicate 
that the departure of the spectral scattering profile from a Gaussian profile 
for an ideal gas nay be as large as 4X at surface atmospheric densities 
(y ‘v 0.3). It follows from Table 4-7 that the effect of Brillouin scattering 
is even leas important at higher altitudes. The relatively small departure 
of the molecular scattering profile from s Gaussian due to Brillouin scattering 
at atmospheric densities will cause an even less significant departure of 
HSRL signal magnitudes from their expecced values. The HSRL system derives 
two spectral integrals from the return signal. The values of these inte- 
grals are insensitive to small uncertainties in the shape of the molecular 
return. This is shown in Figure 4-14 (Section 4-3-g) , where an uncertainty 
in atmospheric temperature of ± results in a measurement error which 

is less than IX. 

Even though the ^impact on HSRL operation is insignificant, the spectral 
intensity departures due to Brillouin scattering can Le accounted for by 
theoretical corrections. Yip (1971) has incorporated kinetic models for 
polyatomic molecules into the Brillouin scattering theory. The results of 
Yip agree with the experimental results of Greytak and Benedek (1966) 
for scattering by polyatomic molecules at atmospheric pressures. 


4- 3-d) Rotational Raman Scattering 

Recent experimental work by Penney et al. (1974) and Fenner et al.(1973) 
has established the absolute cross sections for rotational Raman scattering 
(RRS) from atmospheric gases, particularly for N^, 0^ and CO^. The RRS for 
these simple linear molecules is determined by the allowed transitions 
J J ± 2 where J is the Initial rotational-angular-momentua quantum number. 





/ " 
•A. 

/ 











i The RRS spectral shifts Ao_ T , are approximately 

n 


+ 2 ■ - <« * 6)B o 


4-3, 9a 


- 2 ' + (4J - 2)B = 


4 -3,9b 


where B is the rotational constant for the lowest vibrational level , see 
f O 

1 > \ 

Table 4-8- Values for the lowest possible energy shift * ** B c 

; ! !, 

,are also given in Table 4-8. A comparison of lAo^^l ^ rota Table 4-8 and |60 jj| 

! * } i ’ 

from Table 4-5 showa 'hat RRS from N, , 0, or CO, is not spectrally coincident 
i - [ * i i 

: - ■ i f 

! with the Rayleigh spectral scattering region of interest, and will therefore 

I ! 

not be a concern for the HSKL system. 


Table 4-8 Values of the Rotational Constant for the Lowest Vibrational 
Level for N,,. 0^, and CO^, from Penney et al. (1974) 


B (cm *) 
o 

1.990 

1.438 


i^o-b-al " fryon ) 
11.94 


0.390 


4-3-e) The KSRL Polyetalon Receiver System 

The HSRL receiver system is diagramed in Figure 4-9. To achieve the 
high contrast needed for an operational HSRL, a three Fabry-Perot etalon train 
(FP1, FP2, FP3 in Figure 4-9) is used to suppress background radiation outside 
of its passband of full width at half maximum (FWHM) 60^3 0 ^, the wavenumber 
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of maximum recel ;er transmission. The etalon FP4 is used to achieve very high 

resolution with F »T J; . 6a at a . When a m o , where a is the laser output 

r r o o 

wavenumber maximum, most of those photons scattered by aerosols reach detector 
92, and a large fraction of those photons scattered by molecules (Dqppler 
shifted) reach detector #1. 

As cited in Mack et al. (1963), the instrumental functions* W^ and 
for detectors #1 and 92 are 


W x (o) - w 12 3 <0) ' t " W 4 (o) l + (1 -r t ) 2 * W 4< C >)‘(1~ I^R c) 4-3,10 



w 2 ( o) - w 123 (o). w 4 (o) (l - ~ ,) 8 


4-3,11 


where 


"l23 <0> ' W, 


°* 2 to 123 


a \. (L l Q l do* 


2 S °123 


4-3, 12a 


W A (°) - rM 

4 fin / ; 


o + -Jto 4 


{o V , 

•- rH 




and d ± (o) - (1 + ^ sin 2 (2iTn ' 1^0 cos0)} * 


4-3, 12b 


4-3, 12c 


The notation is as follows: 

A' * etalon coating absorption coefficient 
R’ - etalon coating reflection coefficient 


A The term "instrumental function" denotes the receiver spectral transmission 
function as seen by each detector. 


n' ■ refractive index of etalon gap 
etaion gap spacing 

6 ■ angle of ray path to the optic axis 
i ■ A.(p) ■ the "Airy” function. 

fi °123 acd ^4 are the theoretical resolutions of the PEPSIOS (Poly Etalon Pressure 

1 I 

Scanned Interferometric Optical Spectrometer.) spectrometer the high 

, j * . 

resolution etalon, respectively. The integrals in Eqs. 4-3,12 account for the 

J j ( J 

'decrease in resolution caused by the finite instrumental apertures. Any effects 
due to etalon surface defects have been neglected. 

The maximum angular receiver field of view is given by 48/(1.31R), where 
R - 0/60 Is the effective (actual) resolution of the polyetalon system. The 
solid angle field of view of the PEPSIOS (FP1, FP2, FP3 in Fig. 4-9) Is therefore 


-1 ' )u i23 

ft 123 - 2xR A - 2 * (—Hi) 
r 


4-3,13 


The high resolution etalon must have a smaller field of view to secure operation 
within one free spectral range, such chat 


n - 2rr (— 2. ) 

4 o r 


4-3,14 


The effective receiver field of view can be decreased without a corresponding 
sacrifice in signal luminance by the addition of a receiver telescope of area 
A r « The final "telescopic" receiver solid angle is then defined by the 
constraint that etendue is conserved, such that 

fi r " n !23* a; 4-3,15 


where A a Is the area of the etalon plates, 



I 




4-3~f) Signal Calculation and Solutions 


Let N and H represent the total number of laser photons which are 
backscattered by molecules and aerosols, respectively, and gathered by 
the receiver telescope. Also let N fl represent the total number of solar 
photons which are gathered by the receiver telescope. The quantities N , N 

oi a 

— _ 
and N r are Integrated over the range gated time interval At. The signals 

and S, at detectors #1 and #2, respectively, are linear combinations of 5 , 

* B 

N and N such that 
a s 


S - y. N + y. N + y. N 
1 lm m la a Is a 


4-3, 16a 


S 2 “ Vm + P 2a N a + Vs 


4-3, 16b 


. dN , dN 

y - W q •(— — — * i 

M lm l s e V N do N do' 
m a 


4-3, 17a 


fCO 

i dN ™ . 

m W q • (— — — * — 

2 s e V N do N 

m a 


4-3, 17b 


' - ' : ’ ’ 

W lV ( N~ 

a 


4-3, 17c 


1 dN 

V 2a- M 2V ( IT-d?>* d0 

a 


4-3, 17d 


The operator * denotes "convolution product," such that 


X<0 ) - Y * Z - Y(o) • Z(o-C ) *do 

O J o 


i.'Waf 


f 




^3Vvr 


i ?i^ 


T>N3 


y Is“ W lV d ° 


“2. ‘ j U 2 < >«* da 


4-3, 17f 


The Instrumental functions and are defined In Eqs. 4-3,10 and 4-3,11, 
respectively, and Is the detector quantum efficiency. 


. dN , dN dN 

N~ “do * N - “da and N - ~da re P resent the spectral distributions 
me a 


of photons scattered by molecules and aerosols, respectively, which are 
gathered by the receiver telescope. Since the aerosols are considered to be 
elastic scattering centers, the spectral dispersion due to scattering by 
aerosols 1 b characterized by a Dirac delta function and the spectral distri- 
bution of the laser output is Identically . 

Let the signals measured at detectors ?1 and (2 in the time Interval 

O , A A A * A M 

At be denoted by and S respectively. In general, and 5 * $2 

A A 

due to statistical fluctuations in the photon count* Estimates N and N of 

m a 

. A A 

the quantities N and N can be derived from the measurements S, and S» by 
m a 12 

means of Eqs. 4-3,16 as follows: 


A I At 


N « c. S. + c, S, 
m Jm 1 2m 2 


4 -3, 18a 


N e * c la S l + c 2a S 2 


4-3, 18b 


J 




for 


where 


and 


A • A A 

s. * S. - y n 

1 1 Is s 


h‘ S 2 - u 2. 9 . 


c i„ • + “a . -5 ' 1 


I 


c 2m ‘ - U U‘ 5 


-1 


c la " ~ 


-1 


-1 


c 2« ’ + “to 5 







4 -3, 19a 

fS 

, 

/, 1 

i 

4 -3, 19b 

$ 

* - { 

4-3, 19c 

1 

4-3, 19d 

' 



, 1 

4-3,19e 

-1 


Note that for beat results, we require that the crosstalk terms c 0 and c, 

/m la 

A , 

should be much less than both c^ and c^. *he quantities y^ a N fl and 

A A A 

y N can be estimated ty time averaged valueB of S. and S. during the periods 
Zs s ■*■ ‘ 

between laser firings. 

A A » •• 

The estimates N and N are not equal to the actual values N and N 
m a m a 

A A 

due to two types of error when N and N are derived by means of Eqs. 4-3,18. 

m a 

These are: 

I) error due to a change or uncertainty in an HSRL system parameter. 

Uncertainty in sampled volume temperature (T), uncertainty in laser transmitter 

spectral width (FVfHM ip m 6c m <f.693/b (see Eq. 4-3,5), and uncertainty in 

o o 

transmitter-receiver peak detuning (ip^ “ ° 0 “ <? r ) affect the transmittances 

y <m and y^ a and hence the inversion coefficient* c^ and i * 1,2. 




l i 


- 


-1 




rrv: ... J 
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Type 1 error is calculated below in the following manner: If T, 

A A 

</> o and ^ are perfectly known, then signals and S^Cno statistics) can be 
inverted exactly to original a priori photon inputs and N^. Given change 
in the parameters T, and/or the inverted quantities are calculated 
using the unchanged paraneter values, and the departures of these calculated 

•A A 

values N . N from the a priori values N , N are found. The sensitivity of 

ft D ft Q 

A A 

to uncertainty in T, and is thereby found. 

> , II) error cue to independent statistical fluctuations in the 
photon count at each detector. 

Type II error Is calculated by first finding the total mean number 
of photons incident on each detector. These totals are then given statis- 
tical noise which Is the square root of the total signal, and the corres- 

A A 

ponding maximal departure ofcN., N from the a priori values N_, N_ is 

- < ,y* ft ID ® 


m 


found. This ia sensitivity to number of photons in signal, and this sensi- 
tivity is a function of the ratio $ ■ 5 /N and the a priori solar signal N_. 

am — s 


Type I Error - Uncertainty ia T, f Q and \J> r > , 

, A A 

The following analysis considers deviations of N^, from 5^, N due 

A A 

to an uncertainty in the temperature T. The sensitivities of and N • due to 

uncertainties in the width of the emitted spectral line if^.or the detuning 

of the laser with respect to the receiver ty^.are also given by the following 

analysis by simply replacing T with ^ or . 

o r 

A A 

The departure of N from N is defined by 6N , and the departure of 
mm m 

A A A A 

N from N is defined by 6N . The departures 6N (T) and 6N (T) due to an 
a a a m a 

uncertainty in T can be found by partial differentiation of Eqs. 4-3,18 such 
that 


3c 


6N (T) - S. 
m x 


lm 

3T 


6T + S„ • 


3 ^2m 

ST 


ST 


4-3, 20a 



t v 
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where 6 • N /N . Eqs. 4-3,22 allow us to estimate fractional errors In the 
a m 

measurements due to system parameter uncertainties Independent of the absolute 

values of the signals (S^ - and (S^ “ »2sV* 

The total type 1 error due to uncertainty in all three parameters T, 

4 and 4 can be estimated by 
Y o i . r 


i 

and j 

i 


The detector quantum efficiency q^, and the solar signal N s do not enter into 
either Eqs. 4-3,22 or 4-3,23. 

Type II Error - Photon Statistics 

The number of photons incident on detectors 01 and 02 have mean values 

mm mm A 

and Sj, given by Eqs. 4-3,16. We can expect that the measured signals S. 

/V * 

and reflect statistical fluctuations from these means, such that 


6N 1 ! f SN (T) 2 6N (ip ) 2 6N (4) 2 

— — - (-=__) + ( — — — 2_) + ( — S — L_j 

■ 5 ,• n i S 5 ) 

m | \ m i m m J 


1/2 


4-3, 23a 


6N 1 


A 

[ <— — ) + ( * ---) + ( —! — ■) c 

l «a ' 5 a N a ) 


( '* A A 

6N (T) 2 6N } 2 6N <4 

3 O \ , / 3 O 


1/2 


4-3, 23b 



4 -3,24a 
4-3, 24b 


The effect of these statistical fluctuations on the measurements 
be found by taking the differential of Eqs. 4-3,18 as follows: 


A A 

S and N can 
m a 


« V : - (' 1 ,« V 2 + < c 2 . { V 2 + “ l 2 + ' 2 . “ 2 ^> • ' S “.> 2 4 ' 3 .«• 


“ d <“/ ■ <'1. { V 2 + 4 ('u 4 c 2a ^>-< 5 V 2 4 - 3 > 25b 



■T" 




9b 


The coefficients of the solar terms In Eqs. 4-3,25 are a sum of squeres since 
the statistical fluctuations of solar photons at the two detectors are inde- 
pendent of each other. Substituting Eqs. 4-3,16 and 4-3,24 into Eqs. 4-3,25, 
it follows that 


4-3, 26a 


and 


(6N 11 ) 2 - c? *(y, N + y, N ) + c 2 • (y, N + y 5 ) 

m ' lm 1® m la a 2m 2m m 2a a 


+Ic lm M ls (1 '* 1J l8 >+ C 2m M 2s (l4v 2s )l4S s 


(5N 11 ) 2 - c 2 *(y, N + y J ) + 

a ' la lm m la a 2a 2m m 2a a 


+ i‘u u i.< 1+ V* 


Again letting <J) 'N^, the fractional error in the measurements due to rms 

statistical fluct--ations can be estimated by dividing Eqs. 4-3,26 by N q and N^, 


4-3, 26b 


such that , 


(si u \ 2 

m 


{ <C 1^ »1« + ' 2 « P 2a > + ♦•(c^ Pj. + U 2 „) 


4-3, 27a 


+ f - |c li + 

m 


~ * De 8*(<P) + \ • Deg 2 

N n N 2 mB 

m m 


6N 


II\ 


N / 
a 


■ { ^ c la W la + '2.“2s> * ♦'‘'I'* + W 2 »> 


4-3, 27b 


f - + <: 2 ^ 2 . ( 14 « 2. >1 }- (S .)" 1 


N 


■ r • I * 8 « 2 w *fi" Dee .a 


N 


m 


1 






\ 



96 


- r» — t j— ^ . 1 » » --^v- fr i v *i^ 


The factors Deg m ($), bagJQ) ran be interpreted hs those proportions 

to which the "natural" fractional statistical fluctuations iT^ FT** are 

m * a 

changed (increased) o* the HSRL signal inversion process. The factors 
Oeg^W), Deg fl ($) are both proportional to the Inverse square root of the 
detector quantum efficiency q e> The factors Deg^, Deg^ take into account 
the statistical degradation due to a solar background continuum. 

4-3-g) A Study of Two Example HSRL Systems 

This section presents calculations of HSRL system parameters c. and 

izn 

c , i » 1,2, and estimates of the sensitivity of the solutions to the parameters 

T, - 6cr and ^ ■ a - a . The theory for these calculations was presented 

o o r o r 

in sections 4-3-e) and 4-3-f). The calculations are performed for two HSRL 
systems which differ only in the fourth etalon spectral resolution, see Figure 
4-9. The characteristics of both systems are given in Table 4-9. 

Drawings of the spectral profiles of the laser output and the two receiver 
system passbands as oeen by both detectors are given in Figure 4-10. 

A /\ 

The expected values N and W are calculated from the measurements 

in a 

A A 

S, and S_ by means of Eq6. 4-3,18. The inversion coefficients c. and c. , 
i L lm la 

i ■ 1,2 are functions of T, and ^ only . The inversion coefficient values 
are plotted versus departures of T from 250°K in Figure 4-11, versus departure 
of <1/ from <$a - 1/10 6o„ in Figure 4-12, and versus departures of <(» - l°„ _a r l 

t O O D L u 

from zero in Figure 4-13. All four inversion coefficients are given for both 
example systems. It should be noted from Figure 4-11 that these coefficients 
are relatively insensitive to changes in atmospheric temperature. 

Estimates of Type t error versus uncertainty in T, ip o and are 
given in Figures 4-14 through 4-16. The ordinates of Figures 4-14 through 
4-16 are given in fractional error in measurement aa given by Eqs. 4-3,22. 


Table 4-9 Characteristics of the Two Example HSRL Receiver Systems 


Wavelength 
PEPSIOS FWHM 

High resolution etalon FWHM: 
System one 
| . System two 
Laser output FWHM 

' i r 

t 

I Etalon coating absorption 
’ coefficie'-t 

. , il : 1 ! 

• f 

I Etalcn coating reflection 
coefficient 

1 Detector transmission 
coefficients i 

t 

System one 

\ 

System two 

Solar degradation coefficients^ 
System one 

System two 


4880 X 



e °u3 

- fia D (at 2S0°K) 

■ .086 cm ^ 

(6A n, 20 mX) 

o* 

1 

1/5 fio D = 

.017 

cm * (dA'wfcmX) 

6o. «* 
4 

1/2 6a D - 

.043 

cm (dA'vlOmX) 

i 

o 

D 

«o 

1/10 6a D - 

.009 

cm ^(6X^2 mX) 


A* - .01 
R* - .91 


lm " * 24 

y 2m“ 

.06 

la ” * U 

y 2a m 

.22 

^ - - 17 

U 2m" 

.11 

'la - * 09 

y 2a 

.24 


Deg mp - 5.0 

Deg “ 2.8 
as 

Deg - 6.8 
ms 

Degaa - 


The solar signal transmission coefficients y^ and y are approximated 
by y. and y„ , respectively. 
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Figure 4-17 The ratio i}> ■ N /N as a function of altitude for 
° a m 

the McClatchey et al. (1971) Clear and Hazy model Tropical 

atmospheres. The ratio | is independent of lidar platform 

height h Q . 




Figure 4-18 The HSRL statistical degradation factors Deg (♦) 


and Deg ($) for system one (solid) and system two (dash) as a 
a 


function of the ratio | ■ N /N where q 
. am e 

Devi) « q~ . .... . ... ^ 


Note that 
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1) average slab aerosol extinction coefficient 
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ill) Integrated optical depth 
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v) aerosol volume backscatterlng coefficient 
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vi) average slab aerosol phase function for backscatterlng 
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4-3, 29c 


4-3, 29d 


4-3, 29e 


4-3, 29f 


The albedo for single scattering is approximately 1 at visible wave- 
lengths, but departs significantly from 1 due to 0^ absorption in the UV. 
Typical values for 1 - u)^ as a function of altitude for X ■ 337lX are 
given in Figure 4-19. The data of Figure 4-19 were computed using the 
McClatchey Tropical model atmosphere. Assuming that the albedo for single 
scattering for aerosols and Rayleigh scatterers is a sufficiently well known 
quantity below 10 km, then the above solutions (Eqs. 4-3,29) determine three 

A ’ A 

unknowns (B ,8 >0* 00 ) from two measurements (N , H ). A third 

m a a m a 

measurement is required to solve these equations, and we therefore require 


£ 



/ 


that an independently measured temperature profile T be provided as a function 

of height. As shown In Appendix 4-1. the profile of £ can be determined from 

T in a dry atmosphere with an accucacy —g— < 2 — • The uncertainty in the 

p m 

calculated value for Rayleigh optical depth 


T m <°’*o> 


(X 

o 

e m 


4-3,30 


f 6T m ( 0 ’ x o )l2 * <“? 


f X 

^) 2 -( °6 

i m 

m 


dx) 2 s 4 (^) 2 -[t (0.x )] 2 

I m o 


vpnj 

where It has been assumed that ( — g — ) is a constant independent of height. 


L molvcutor olb#do for 
ttn£f« tcoflCnnQ 


I McCLATCHEY TROPICAL ATMOSPHERE 

| 40 L v 3371 A 


Figure 4-19 Typical values for the molecular albedo for 
single scattering 1 - w as a function of altitude.from 
the McClatchey et ai. (1971) Tropical model atmosphere 
at 3371 X. 
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The errors are as follows: 

1) average slab aerosol extinction coefficient (.from 4-3, St) 

i«b^ )■ ■ -- 1 - ■ , (m. (o.vi* ♦ nt. w-V 1 ’} 

— Sj ) 

j . i 

iij aerosol optical depth (from 4-3, 29b) 
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4-3, 31b 


ill) integrated optical depth (from 4-3, 29c) 
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iv) scattering ratio (from 4-3, 296) 
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v./ aerosol volume backscatterlng coefficient (from 4 -3, 29e) 
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vi) average slab aerosol phase function for backscatterlng (from 4-3,29f) 
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Substituting Eqs. 4-3,30, 4-3, 31b and 4-3, 31c into 4-3, 31a, and letting 
60 * T 

— - 2(—=r) and (x. - x.) - u, then Eq. 4-3, 31a can be written 
p x i j 

m - 

. , f 6 y .» « 6u> « 2 dTll** 

« r 1J - (_= > 2 ♦ <-f ) 2 + <.-f ) 2 + w + ‘' 0 2 < o .« 0 )! • < t > } 

a uV2 [ ’m m m J 


6S , 6N 1 , 6N 11 

where we substitute ( — — ) - ( ) + ( ) from section 4-3-f). The 

N N N 

' - ® m ■ n 6Y m 6w n 

factor T (o,x ) can be neglected since 4t (o,x ) « 1. Letting — — - “ — — 
mo mo i w 

A mm 

6n 6n 6y 6t 

then ^ is a function of u, — — (or — — ) , — — and , such that 

4 S N y m T 

m m 


v a TT 

~ 6 N A „ 6 N Li 


( dF ^) 2 - -K K --) 2 + ( — ) 2 + 2 (- r ) 2 + 4 (^) 2 ] 

a r- r. T 


4-3, 32b 


2u . N 


The HSRL instrumental error , given by Eq. 4-3, 23a, 1 b due 

- N >■ - 

O _ 

to the uncertainry in the parameters T, <j> o and Partial errors due to 

uncertainty in T, and ^ are given for the example HSRL systems in 
Figures 4-14 through 4-16. For example, these figures show the following 


errors for Systems one and two (4880X) : 


6N (T) 


- 0.10% (one). 


0.15% (two) 


for <$T - + 5*K 


«»(*„) 
m o 


2.0% (one). 


-0.5% (two) 

~ "1 / t 


for 6^ -5 x 10* cm” 1 (= 0.1 mX) 
o 


6N(iJ») 
m r 


-0.5% (one), 


- 0.1% (two) 


for 6* -5 x 10' 4 aa _1 (= 0.1 mX) 


l 


i - . 
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The HSRL calibration constant can be determined accurately, such uhat 


6y 


* 0.5Z . 


'm 


the uncertainty in atmospheric temperature is 
6T 


= 1.0% 


j_4 

The quantity 6$ a is given for various values of u and — — in Table 4-10, 


6N 


II 


given the above values for 


«V 

K ’ Y * 


m 


and 



Table 4-10 Typical Values for 60 

(km ’) as a 

Function of 



6N 
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u * 100 m 
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500 m 

1 km 

i% 


.235 km" 1 

.117 km -1 

.047 km _1 

.023 km" 1 
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.308 
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.062 
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10 % 


,742 
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4-3,33 


Error due to statistical fluctuations is a function of the ratio $ ■ N /N , 

8 OB 

the number of photons which are backscattered by molecules H , and the number 

ID 

of solar photons N fl . The ratio 4> can be determined as a function of height 

from Figure 4-17, the total number of collected photons which are scattered 

by molecules can be determined from either Figrres 4-3 or 4-4, and the solar 

signal can be determined from Table 4-5. The statistical degradation factor 

Deg(<|>) can be determined from Figure 4-18. The factor Deg is given in 
, m ns 

Table 4-9. j The error due to statistical fluctuations is then 


lotting N » — — • At , the signal integration time required to obtain a 
s m at m 

specified statistical precision in the molecular return is 


At “ 
m 


Deg 2 (<J>) + (N /N ) • Deg 2 
m s m "ms 

K'V 


4-3,34 


The type II or statistical error 
6(T dJ from Lq. 4 -3,34. 


can be determined as a function of 


- 3 

\ 

t • “ • 

> : 


The precision to which 60 must be determined lr. a 1 b vertical Interval 

a 

in order to obtain an accuracy of ♦ 0.2 °K day ^ and ± 0.5 9 K day ^ in the total 

daily heating over land was estimated in Tables 3-1 and 3-2, respectively. These 

results are used to obtain estimates of the HSRL signal integration time for the < 

model HSRL systems one and two (Table 4-9) located at h ■ 185 km (Shuttle caae) 

° l 

and h ■ 10 ka (aircraft case). The signal integration times are given as a function! 
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of height in Table 4-11 for the cases of a McClatchey Tropical Clear (BF ■ 1) 
and Hazy (BF ■ 5) model atmosphere. The signal integration time-. &t ffi , ia 
most sensitive to the statistical degradation factor Deg^ (<p ) . Thf HSRL model 
system one is defined to a greater spectral resolution than system two, and it 
has a correspondingly lower value for Deg n (4>). System one can therefore oper- 
ate at lower laser power levels to attain a specified statistical precision in 
the molecular return signal. The calculations of Table 4-11 were performed 
using the worst case solar background as given in Table 4-5. The solar back- 
ground was reduced an order of magnitude for the Shuttle case to characterize 

HSRL operation within a Fraunhofer line. The values for At calculated with- 

IB 

out a solar signal are also given in Table 4-1! . The length of this sample 
path on the Earth's surface, assuming that the Shuttle orbit is circular and 
that the aircraft speed is 800 km/hr, is also given in Table 4-11. 

4-4 Lldar Methods to Measure Aerosol Profiles 
4-4-a) System I: HiRh Spectral Resolution Lldar 

The High Spectral Resolution Lldar (or HSRL) concept was presented in 
Section 4-3. An envisioned HSRL system would employ an N_ laser to pump a dye 
cell, which is tuned to yield a narrow spectral line by means of Fabry-Ferot 
etalons. A high resolution receiver which employs Fabry-Perot etalons is also 
required. Independently derived temperature profiles are required to provide 
vertical profiles of the Rayleigh scattering cross section. Eqs. 4-3,32 
suggest that temperature profiles derived from the inversion of data from 
passive satellite borne sensors will suffice for this purpose. 

Source : Nj superadiant (337lX) laser pumping a dye (4000-5000 X) in a 

laser-etalon-amplif ier configuration with an on-axis etalon. 

Spectral width of 0.002 X. Energy per pulse - 50 pJ. Repetition 
rate “ 60 Hz . Pulse duration ■ 7na . Energy conversion efficiency 
■ (0.1% N 2 conversion) x (30% N^-dye conversion) - 0.03% . 
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Table 4-11 Signal Integration Tine and Sample Path Length for Model HSRL Syates ONE Located at Shuttla and 
Aircraft Altitudes. ^ 



t 0.2*K day* 1 Accuracy* 2 * 


t 0. 5*K day 

-1 . (2) 
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\ O — 
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J 

2 

47.0 

355. 
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140. 

.372 
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3 
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.235 
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.635 

141. 

j 
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.781 
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.445 
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J 

5 
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Table 4-11 (coot.) 

(1) lb* model HSRL systems are defined 1 a Table 4-9 (l * 4880 X). Additional nodal 
parameters are e« followat 

9, • 0.25 { 

l f • 1.0 n 2 (Shuttle) 

• 0.1 n 2 (Aircraft) 
t 0.2*1 day” 1 accuracy! 

C - 10 bJ sec” 1 (Syaten one) 


! w 


• 20 aj aac 

• 1.0 ion 
-1 


-l 


(Syacaa two) 


± 0.5*K day ~ accuracy: 

E “ 3 «j aec 1 (Syacaa one) 
o 


- 6 aj ecc 
• 0.75 km 


(Syacaa two) 

/ 


# 60 Ha repltltion rata' 
€ 60 Ha repltltion rata 


g 60 Hi repltltion rata 
? 60 Hi repltltion rata 


The veret caao aolar algnal fro* Table 4-5 la Included. 

: H : ! 

Values for | A3 | are taken froa Tables 3-1 and 3-2. and correspond to an uncertainty 
*1 .1 

of 1 0.2*1 day and t 0.1*1 day In the total dally heating over land. The model 

atmosphere parameter BF - 1 end B? - 5 correspond to the KcClatchey Clear and Bray 

nodal ataoepheras, respectively. The values for I A0 I and p - H /S are as follow: 

• an 




Altitude 


KcClatchey Tropical Clear 

A I AO I * Ilf 1^ 


KcClatchey Tropical Harr 


(km) 


(km-1) 

(km-1) 


(to" 1 ) 

(to" 1 ) 

C 

2.4 

.08 

1.1 

10. 

.13 

.25 

•- * 1 

1.3 

.06 

.08 

4.4 

.43 

4.3 

2 

.62 

.03 

.11 

1.8 

.09 

.23 


.31 

.06 

.10 

.73 

.13 

.33 

4 

.18 

.04 

.13 

.29 

.03 

.19 

5 

.13 

.03 

.30 

.14 

.11 

.43 

* 10.2*1 day' 1 

accuracy 






b 10.5*1 day' 1 

accuracy 
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( 3 ) 


For the Shuttle case, the sample path lergth along the Earth's surface Is As - 

where the angular velocity for a circular orbit Is 8 ■ — — 

a3 


and a - 


*e 


+ h 


- mean Earth radius - 6.4 : 10 b 
C ” gravitational constant - 6.67 x 10 

• nass of tha Earth ” 5.98 x 10 
* -1 


-20 


24 


. 3 -2. -1 

ka aac kg 


Cor h - 185 ka. 4 ■ 7.56 ka sec 
o 

(4) Tha aircraft velocity over tha Earth's surface la assuaed to ba 800 ka hr” 1 . 

(3) Solar signal values are derived fron Table 4-5 (1 - 4880 X). Solar signal values 
for tha Shuttle case vara reduced an order of sagnltuda to characterise HSU. oper- 


ation within a Fraunhofer line. 
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Receiver ; Four etalon system for high resolution (0.002 X) notch and spike 
and high rejection of scattered solar continuum. Transmittance 
of instrument 3 35%. Quantum efficiency of photomultipliers 3 30X. 

Overall Efficiency - .0031. 

In order to achieve stable wavelength outputs as narrow as a few milli- 
angstroms at high power It appears necessary to use low feedback in the 
laser system In order to avoid being locked into laser cavity modes. The 
wavelength oi the cavity modes depends critically on Inhomogeneities in the 
dye cell and the separation of the laser cavity mirrors. The separation 
needed between laser cavity mirrors is generally so large that a linevidth 
much narrower than needed is produced. 

Waller8tein and Hansch (1975) recently re, rted that an Nj pumped dye 
laser-amplifier system is capable of linewidths as narrow as 0.006 cm * 

(1.5 mX at 4800 X) with pulse energies of 500 yj when pumped by a 10 m j , 

10 nsec pulse from the laser pump. Projecting this conversion efficiency 
to the case of a 0,3 watt laser (average power for 10 nsec pulse lengths 
X a 60 Hz repetition rate), it appears that high resolution (1.5 mX) output 
from the dye laser-amplifier may be feasible with an average power of 15 mw. 
This would be adequate for the proposed HSRL system. 

The envisioned laser system is configured as follows: The original 

design of Hansch (1972) for en pumped dye laser is employed. This design 
consists of a dye cell enclosed in a cavity which includes an uncoated flat 

i 

end window, the laser pumped circulating dye cell, a beam-expanding tele- 
scope, an intracavity etalon and a high blaze-angle echello grating end 
reflector. The output of the dye laser is coupled through a wide-spaced 
confocal Fabry-Perot etalon into two amplifier dye cells in succession. Each 


i 


~~ 1 ” 
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dye cell does not require a resonant cavity rnd each Is pumped by a fraction 
of the same laser pulse which initiated the dye laser operation. The 
output of this system has a spectral width determined largely by the width 
of the transmission peak of the external wide-spaced etalon, and by the 
effects of Heisenberg broadening due to the extremely short duration of 
pump excitation. An attractive feature of this system is the absence of 
cavity mod. »s other than those of the external resolving etaloa, and precise 
tuning of the laser-amplifier output to the high resolution PEPSIOS receiver 
appears feasible. 

\ . . 

Preliminary work has been carried out at the University of Wisconsin 

Physics Department to investigate the possibility of using such a laser- 
amplifier system for the HSRL project. Several dye laser-amplifier con- 
figurations have been tested. The most promising configuration employs an 
Nj pumped dye laser of essentially the original Hansch design, with the 
output coupled through an air-spaced plane Fabry-Perot etalon to a single 
dye amplifier cell. The output of this dye laser amplifier configuration 
is shown in Figure 4-20. The resolving (external) etalon used an etalon 
spacing of 9 mm with a coating finesse of 30." An amplified output of 
3.5 mwatt (average power with a 10 nsec pulse length at 60 Hr operation) 
was obtained using a 0.12 watt (average power) Nj laeer pump. The spectral 
profile of the output closely follows the transmission profile of the 9 mo 
spaced resolving etalon which naa a measured width of 5 m? at 4800 £. We 
expect that future work with a resolving etalon of wider spacing (higher 
spectral resolution) will allow a reduction of the spectral width to the 
requirements of the HSRL system. 




i-i 


A 


u ctaifat 



Figure 4-20. High resolution interferometer scan of N 2 pumped 
dye. laser output, showing spectral width of 5 mX (FVfHH) at 
4880 X. Average power of the amplified output was 3,5 mwatt 
with a 10 nsec pulse length at a 60 Hz repitition rate. 

We have not had an opportunity to investigate flashlamp pumped dye sys- 
tems to determine whether our requirements could be achieved with them. 

Their generally much longer pulses (100ns or more compared to 'v 3 ns) may 
make it difficult tc attain efficient operation while avoiding ultra-narrow 
but unstable cavity modes. Generally, the overall efficiency of conventional 
flashlamp-pumped dye-systems is higher than for ^ laser-pumped systems. 

1 

1 

Some considerations regarding this technique: ./ 

Pro: 

1) The aerosol extinction coefficient 8 is uniquely determined at a parti- 

a 

cular wavelength. 













2) The instrumental requirements have been demonutrated in the laboratory. 

3) The science of laser remote sensing of the atmosphere by means of ground 
and aircraft based lidars is promoted by research in this area. 

4) Solar heating rate computations require information about over a broad 
wavelength range. If raultispectral imagery with high spatial resolution is 
available, then the technique described in Appendix 1V-2 can be used to derive 
the total optical thickness of the aerosols as a function of wavelength. If 
it Is assumed that the aerosol size and composition are independent of alti- 
tude then B can be inferred as a function of wavelength for G.3um <X< 2.5um. 

a 

Passive polarization measurements such as those proposed for the nimbus G-ERB 
experiment by Stowe aud Hilleary (1975) may also provide information on 
aerosol sizes. It should be noted, however, that the analysis of such data 
may be ambiguous (see the remarks by Kuriyan in the reference cited above). 


1) Infrared heating rate computations require prior knowledge regarding 
aerosol composition and size distributions before the 0 a values derived at 
solar wavelengths can be introduced into a GCM. This problem is common 

to all of the lidar techniques which are presented in this section. 

2) Return optical power levels are marginal at Shuttle Tanges. However, 

high resolution laser technology is expected to advance over the coming years. 

3) Receiver optics is sophisticated, and may require frequent maintenance. 
Engineering developments should assure reliability by the time that the 
Shuttle mission will occur . 

These considerations render the KSRL technique most attractive for ground and 
aircraft applications. The HSRL technique is somewhat marginal for Shuttle 
applications, but the present limitations may become less compelling in the 


near future 





4- k-b ) System II; Multi-Wavelength, Low Spectral Resolution Lidar 

This system employs an array of three or four flash lamp pumped dye 

lasers, fired sequentially at a high repetition rate. It may alternatively 

employ an laser to pump four dye cells. Photomultipliers measure returns 

i ■ 

at each wavelength Xj(J “ 1,...,4). Each wavelength band is defined within 

6X - lX resolution. ’ - ' ‘ * 

i ■ ! 

• . *■ i *' ' ' j 

Source: Multiple wavelengths by sequential firing of flash lamp pumped 

\ dye cells. Alternatively, multiple wavelengths by superposition 

'f, ‘ of output from multiple dye cavities pumped by a common N z laser. 

, . , Spectral width of each line * lX. Energy per pulse (all wave- 
lengths) «• 50yJ . Separation of lines - 1000 A in the range from 

4000 to 7000 X. Repetition rate *■ 1 Hz (flash lamp), «■ 60 Hz(N z ). 
Pulse duration - 20ns (flash lamp), » 7ns (N). Energy conversion 
efficiency - 0.3% (flash lanp) , = 0.03% (N z ) 

Receiver : One interference filter (1 X resolution) per wavelength. Filter 

transmittance = 50%. Average quantum efficiency of photomulti- 
' pliers=30%. 

Overall Efficiency : - 0.01% for pumped system, 0.1% for flash lamp 

pumped system. 


Theory: 


This scheme measures - 

rx 

h.s ‘ifr- < "?.l \,s + ! 6 ».i < s .,j + 

° ,J o 
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Four such measurements are taken at X^ (1 < j < 4). The Rayleigh extinction 
profile is independently measured (e.g. from radiometric inversions). The 
aerosol backscattering phase functions and the extinction coefficient are 
assumed to be of the form: 


IP* 11 } - AX? 

a.j J 


4-4, 2a 


8 - CX? 

J 


4-4, 2b 




and 




" p xl j. 

$h 




Thus the method yields four numbers for each layer, viz A,B,C,D. If the 
aerosol size distribution and composition are slowly varying functions of 
altitude, some of these numbers may be independent of altitude and the equa- 
tions will be overdetermined. 


1) The technology is relatively unsophisticated. Flash lamp pumped dye 
lasers have been flown on aircraft. 

2) The power utilization efficiency of flash lamp pumped systems is approxi- 
mately .J.t, which is convenient for remote stations with restricted power 
sources. 

3) An estimate of aerosol size can be derived directly. 


1) Prior knowledge of the optical properties of aerosols Is required to 

constrain the choice of solutions to power laws; i.e. the representation 

0 . • CA® assumes that the aerosol size distribution obeys a power law. 

0 1 J J 

2) Prior knowledge regarding aerosol composition is still required to Incor- 
porate results Into GCM. -j-'V 

3) The lifetime of the flash lamps or optical components subject to high 
peak powers may be restrictive. 

4) The simultaneous operation of four laser systems may be operationally 
difficult. 

5) The science of remote sensing of the atmosphere by means of lasers wilj 
not be significantly advanced. 

4-4-c) System III: Angular Scanning System 


* 


A lidar system is considered which alternately fires toward nadir and 
at an angle ahead of the vehicle which transports the lid-r. The same volume 
of atmosphere is thus probed at two angles at slightly different times. 


/ .3 ■ 



/ 




m 
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DeLuiai et al. (1975) have explored a ground based Ildar system which utilizes 
a combination of the multiple wavelength end multiple > ngle approaches. 

1 " V 

Theory: 

Given a Ildar system which measures N Q «* N + N from the same atmo- 

0 m a 

spheric volume at angles 0^ and 6^, the lidar equations to be solved are 


, IP (*) 

N(0. ) - — (^r 0 to + -rz — 0 to ) • exp [-2 sec0. (0 + 0 )dx'J 4-*, 3 

it 2 q in to fl fl k. m a 

x sec0. 
ok 

o 

E Q A r n r u 

for k ■ 1,2 and Yg ■ — — . Assuming th3t the range interval u is 

independent of pointing angle, and that the atmosphere is horizontally 
homogeneous on the average in the region of interest, the various optical 
properties (analogous to Eqs. 4 - 3 , 2a) can be solved for as follows: 


i) average slab aerosol extinction coefficient 


jpij _ T a (0,x i ) ~ 

P a (x £ - Xj) 


4-4, 4a 


il) aerosol optical depth 


T (0,x ) » T(0,x ) - 0 dx 

a o o m 


4-4 ,4b 


..*3- 


iii) integrated optical depth 

^ sec0^ • N(0^) 

t(0,x q ) 2 [sec 02 - sec0^] secf^”"*"" IT(0^) 


4-4, 4c 


iv) scattering ratio 


N(0) x 2 sec0 
o 

V 8? 8 «“m 


• exp 2secG 


• T (°* X o> ] 


4-4, 4d 


v) aerosol volume backscattering coefficient 


B (tr) - 

-4 - (s - i) — B u) 

/in 8ir m m 


4-4 ,4e 
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vl) average slab phase function for backscatterlng 




-Ij 


4ti 


B a <*> —u _! 

— - [oj F^ J i A 
4ir 1 a ^a J 


4-4, 4f 


It is very Important In the error analysis that the quantity <t - jsecB^-eecOjJ 
should be as large as possible. However, the equations are physically unreal- 
istic as 6. or 0_ approaches n/2 radians. A good selection would be 0, • zero 

i — * i — ? 1 

1 1 ! a .. . 

radians (vertical) and 0„ - it/3 radians so that $ • 2. 

!■ , | 

* j * 

The errors are as follows: 

i) average slab aerosol extinction coefficient (from 4-4, 4a) 


nr 11 ! 2 

a 




2 j(6T a (0,x 1 

- X.) ( 


)] 2 + [ fi ^^. Xj )] 2 


} 


4-4, 5a 


li) aerosol optical depth (from 4-4 ,4b) 

[6t a (0,x o )] 2 - [6 t(0,x o )] 2 + [6 t b (0,x o )J 2 


4-4, 5b 


ill) Integrated optical depth (from 4-4 ,4c) 


[6T(0,x o )r - 


4(sec0 2 _ secQj^) k-1 [ 


6N(0 k ) 


W 


4-4 ,5c 


iv) scattering ratio (from 4-4, 4d) 


63 


(^f) 2 - [~^] 2 (|^p) 2 + exp(4sec6 • t(0,x n )J* [6 t(0,x / i )] 2 4-4, 5d 


N(0) 


v) aerosol volume backscattering coefficient (from 4-4 ,4e) 


B (it) 2 

f — - — r 
1 B a (^) ] 


<f > 2 * <> 2 * £> 2 

m m 


4-4, 5e 
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vi) average slab aerosol phase function for backacattering (from 4-4 ,4f) 


6 IP _ (it ) ij 6B (it) , 6F lj , &•» , 

( -- - ) “ < i — oPr ) + ( > + t"7T' 

IP fir) 1J B a (lt) ^-lj w a 

a 3 a 


4-4, 5f 


Pro: 


1) The aerosol optical properties are well defined at a particular wavelength. 

2) The laser and detection systems could be unsophisticated, l.e. low resolu- 

tion flash lamp pumped dye lasers and low resolution receiver systems could 
be employed. >" : 

• " __ * __ - *■ ' f • 

Con : , f ^ ‘ • 

1) Figuring aim. diameter receiver mirror to provide good optics at *v ± 30* 
off axis appears to be expensive. The implementation of this scheme on a 
Shuttle appears to be more difficult than the preceding schemes. A smaller 
rotating mirror could be mounted in a pod outside an aircraft, but the cost 
of building an F.A.A. approved pod would be high. A multiple transmitter 
system is also a high cost item. 

2) Prior knowledge regarding aerosol size and composition is required to 
incorporate the 6 profiles derived by this technique into a G.C.H. 

») Although considerable technical effort would be required for mirror 
engineering, the science of remote sensing by lidar from a Shuttle would 
be negligibly advanced. 
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4-4-d) System IV: Broad Band High Resolution System 

This scheme Is similar in principle to System I (HSRL) except that all 
of the light from an N_ laser is utilized directly. The broad spectral out- 
put of the laser is spread in angle as a function of wavelength so that high 
spectral resolution is obtained as a function of angle. The receiver system 


. J 

* i 

' | 
i 


is matched to the transmitter so that radiation which is received at a par- 
ticular angle must arrive with the appropriate corresponding wavelength, see 
Appendix IV-3 . 

Source: N laser-etalon-amplif ie'r combination with off-axis etalon to give 

a wavelength-encoded angular distribution of laser light with a 
recolution of 0 .001 X at 3371 X. Energy per pulse ■ .003J . Repe- 
tition rate « 60 Hz . Pulse duration ■ 7na . Enerry conversion 
efficiency 3 0.1% . 

Receiver: Off-axis three etalon system for high resolution (O.OOlX) notch 

and spike and moderate rejection of scattered solar continum. 
Transmittance of instrument “ 3%. Quantum efficiency of photo- 
nultiplier ■ 40%. 

Overall efficiency 3 0.001%. 

Limitation : Only useful over short ranges (ground or airplane based) due to 

possible multiple scattering across wavelength encoded zones in 
the field of view, which would effectively degrade the effective 
spectral resolution of the system. 


Pro: 

1) The derived quantities are well defined at a particular wavelength. 

2) The techniques have been demonstrated in the laboratory and the system is 
reasonably straightforward to implement, (see Appendix IV-3a.) 

3) The power utilization efficiency lies between methods #1 and #11. 

4) The science of laser remote sensing from the ground and aircraft can be 


advanced . 











123 


References 


ANSI, (1973), Safe Use of Lasers . American National Standards Inatltute, 
Institute, 7-136.1, New York. 

Barmore, F.E., (1975), The Filling-In of Fraunhofer Lines in the Day Sky, 
J. Atm. Sci., 32 , pp. 1489-1493. 


Cummins, H.Z., N. Knoble, and Y. Yeb, (1964), "Observation of Diffusion Broad- 
ening of Rayleigh Scattered Light," Phya. Rev. Letters, 12, 150-153. 


Dabberdt, W.F., and W.B. Johnson, (1973). "Atmospheric Effects Upon Laser Eye 
Safety — Part II, Final Report Stanford Research Institute, Project 7472. 

DeLuisi, J.J., B.G. Schuster and R.K. Sato (1975); Separation of Dust and 

Molecular Scattering Contributions to the Lidar Observation: A Method. 

App. Opt. 14, 1917-1923. ... 

Fenner, W.R. , H.A. Myatt, J.M. Kellam, and 3.P.S. Porto, (1973), "Ranan Cross 
Sections of Some Simple Gases," J. Opt. Soc. am., 63^, 74-77. 

Fiocco, G., and J.B. DeWolf, (1968), "Frequency Spectrum of Laser Echoes From 
Atmospheric Constituents and Determination of the Aerosol Content of 
Air, J. Atm. Sci., 25, 488-496. 


Fiocco, G. , G. Benedetti-Miche] angeJi, K. Maischberger , and E. Madonna, (1971), 
"Measurement of Temperature and Aerosol to Molecule Ratio in the Tropo- 
sphere by Optical Radar," Nature, Phys. Sci., Letters, 229 , 78-79. 

Grainger, J.F., and J. Ring, (1962), "Anomalous Fraunhofer Line Profiles," 
Nature, 193 , 762. 

Hansch, T.W., (1972). Repetitively Pulsed Dye Laser for High Resolution 
Spectroscopy, App. Opt. 11 , 895-898. 

Kitt, Peak Solar Atlas (1975), in press. 


Mack, J.E., D.P. McNutt, F.L. Roesler, and R. Chabbal, (1963), "The PEPSIOS 
Purely Interferometric High-Resolution Scanning Spectrometer. I. The 
Pilot Model," App. Opt., 2, 873-885. 


McClatchey, R.A., R.W. Fenn, J.E.A. Selby, F.E. Volz, and J.B. Garing, (1971), 
"Optical Properties of the Atmosphere (revised)" AFCRL-71-0279. 




R 




i 

3 


Moses, R.A., Ed., (1970), Adlers Physiology of the Eye; Clinical Application . 
5th Edition. St. Louis, C.V. Moaby Co., Chapter 11. 

Hunter, D.M., (1970), "Surface Albedo and the Filling-in of Fraunhofer Lines 
in the Day," Astrophys J., 159 , 1107-1110. 

Noxon, J.F., and R. Goody, (1965), Atmos, and Oceanic Phys., 163. 






124 


Penney, C.M., R.L. St. Peters, and K. Lapp, "Absolute Rotational Raman Cross 
Sections for N 2 , 0 2 , and C0 2 ,” J. Opt. Soc. An., 64, 212-216. 

Pressley, R.J., ed., (1971), CRC Handbook of Lasers , The Chemical Rubber Co., 
Section 1, "Ocular Hazards" by A.M. Clark, pp. 3-10. 

Stowe, L.L., Hilleary, D.T. (1975) Proposal: Adaptation ->f ERB Scanning 

Channels for Polarization Measurements. 

Haller stein, R. and Hansch, T.W., (1975), Powerful Dye Laser Oacillator- 
Amplifier System for High Resolution Spectroscopy, Optics. Comm. 

(to appear July 1975) 

Wang, W., and G.A. Domoto, (1974), "The Radiative Effects of Aerosols in 
the Earth's Atmosphere , " J. Appl. Meteor., 13, pp. 521-534. 

Yip, S., (1971), "Rayleigh Scattering in Dilute C-ases," J. Acous. Soc. of 
Amer., 49, (Part 3), pp. 941-949. 

Yip, S., and M. Nelkin, (1964), "Application of a Kinetic Model to Time- 

Dependent Density Correlations in Fluids" Phys. Rev., 135 , pp.A1241-A1247. 




V.' 


125 


Appendix IV-1 

On the Derivation of Atmospheric Penalty from a Temperature Profile i n a 
Dry Atmosphere . 

Using the equation of state for an Ideal gaa and the hydrostatic approxi- 
mation, the derivative of pressure P with respect to height z is 

/• • 


dP 

ds 


- ET - £T 4*T ♦ Kp P 


ds 


ds 


IV-1,1 


where R j* 2.87 x 10 6 erg gm" 1 °K 1 is the gas constant for dry air and g is 
, the acceleration of gravity. Eq. 1 has the solution 

J * 

! P(*) ■ P(h) exp l-J* f (z)dz] IV-1, 2 

' 1 

for 

f(z) *» T' 1 • g + gR _1 *T -1 

dz IV-1, 3 

where p(h) and T(h) are the atmospheric density and temperature at height h, 

respectively. The error in this calculation is 


V h S(h) ; 


2 + ‘^fuK ♦ i****' 




1 (z)dz) J 2 


IV-1 ,4 


+ ( «IM) 2 

1 T(h) 1 

Breaking up the las; term of Eq. IV-1, 4 into N layers of width u, and letting 
~ » be constant throughout the atmosphere, then 
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Appendix IV-2 : Determination of the Dependence of the Optical Thickness 

of Aeroeolo on Wavelength 

The total optical thickness of aerosol layers can be determined as a 
function of wavelength by means of passive imagery of the earth's surface. 

This appendix will outline a possible technique that may be employed for 
this purpose. 

IV-2. 1 Introduction 

The contrast of images of the earth's surface obtained with spectro- 
photometric instrumentation carried on board aircraft and earth orbiting 
satellites is frequently degraded by atmospheric haze. This problem has 
been comprehensively reviewed by Duntley et al. (1964). The solution of 
radiative transfer problems has advanced since Duntley et al. (1957) 
identified the factors which affect image contrast degradation so that 
acceptable models of contrast degradation can now be developed. 

A method is presented to compute the degradation of contrast as a 
function of the following parameters: 

1. Optical thickness of the haze 

2. Scattering characteristics of the haze particles 

3. Solar zenith angle and viewing zenith and azimuth angles 

4. Mean surface albedo 

5. Contrast of surface features 

The following analysis may be used to determine the optical thickness 
of the haze by measuring the radiances reflected from small objects. 

This problem has been investigated by Breitling and PilipowskyJ (1970), 
Van de l!ulst (1971) and Kondratyev et al. (1973). These studies considered 
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surface features which are horizontally seal-infinite i.e. the mean free path of] 
photons is small compared to the scale of surface features. Sucn schemes are ap-/ 
plicable to surfaces which have features on a scale which exceeds M. km. 

However, much imaging over land masses concerns itself with features 
which are on a scale which is smaller than M. ka. The present analysis 
addresses Itself to the computation of contrast of such features viewed 
through hazes. 

1V-2.2 Analysis 

Consider a plane parallel hazy atmosphere of optical thickness, il- 
luminated from above by a parallel solar irradiance,F o ,per unit area perpen- 
dicular to the incident beam. The radiances diffusely reflected and trans- 
mitted are I and I respectively and they can be expressed in ster ^ for unit 
R T 

incident irradiance in a given spectral pass band. The reflection and trans- 
mission functions, R and T are defined in terms of the radiances in the 
manner described by Irvine (1968)? 


- 


I„(A,,A, y,U„,<» - 


R f 


l (1 + 6 M ) cos 

m-o B »° “ ° 


IV-2 ,1a 


I x (A f .A.M.y o .<}») -- 2 -P o f o 


m,0 


J 0 (1 + 6 n.0 )T o {M ’ P o )c08 

if m -0 
if in i 1 0 


IV-2, lb 


Where 8 Q m cos 


-1 


-1 


U D and 9 “ cos p are the solar and viewing zenith angles re- 


spectively, and $ is the azimuth angle measured from the solar azimuth shown in 
Pig. TV-2,1. 


*The adding procedure is described in detail by van de Hulst (1962), A New Look 
at Multiple Scattering Mimeographed Report N.A.S.A., Goddard Institute of Space 
Studies, New York. Unfortunately this report was not widely circulated. 
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The surfaces are assumed to Lambertian so that radiances reflected froa 

the surface are azimuthally symmetric. Specular reflectors would contribute to 

azimuthally dependent terms, R (p , p ). The underlying surface is assumed to 

m o 

be characterized by an albedo A. for the small features and by a mean albedo A 

of the surface, averaged over a horizontal scale exceeding the order *\* 1 ka. 

The reflection and transmission functions of the haze alone are designated 

R* and T*; they can be found by the doubling method which is described by 
mm 

Hansen <1969), for example. The reflection function R la found by adding 
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Ctoa pairs of surface features measured at a common zenith angles la inde- 
pendent of azimuth angle, viz 


* 1 W A> \i,u o ,0-\{k 2 ,k, y,y o ,4»>] 


= ^(A,T,y,y o ) 


IV-2,8 


l W A 1 ' V 

It an observer has prior knowledge regarding the albedos of recognizable 
features A^, and A and if he measures M q and y, the difference In measured 

radiances can be used to determine the optical thickness, T, of the atmospheric 
haze trom the ip function. 

The optical thickness of a haze can be determined if no prior knowledge 

regarding any surface features is available. The difference, 6l ,, between 

R,t 

radiances reflected from surface features which may be both shaded and unshaded 
from direct sunlight by a large object, can be computed from IV-2,6: 


61 


-t(— + — ) 


R»1 


—t(— ; i) 

, - 5 Vo^ » "o' 

61 - — — e 


IV- 2, 9 


R,2 11 

These quantities are proportional to the solar irradiance which is reflected from 
the surface feature and which passes through the haze undeflected by scattering. 
Combining eqs. (8) and (9) yieldsan expression which is a function of optical 
depth, the mean surface albedo, and the viewing and solar zenith angles only: 

[I (A , A, u, A * M ’ t( u + if - * - 

[ ft- ~ 6I R f — " e U P o ^(A,T,V,y o ) IV-2,10 

R , 1 K,2 

The optical thickness of the haze can be obtained from this expression because 
y and p Q can be measured and A can be estimated. 

IV-2,3 Numerical Examples : 

Various mixtures of black earth and deciduous trees viewed at wave length 
X « 0.8 ym are modelled to illustrate the degradation of contrast produced by 
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h*s«. lhis wavelength la selected for the exaople because Rayleigh scatter las 
and gaseous absorption nay be conveniently neglected, however the present I 

method la applicable to other wavelengths if due attention ie devoted to j/ 

i 

these considerations. 1- 

! 

The hare is assumed to be represented by a water haze K defined by (. 

S'- . 

Deirmendjian (1969). This model assumes that the haze particle size distri- 

V ' 

bution is represented by: 

1/2 

n(r) - 5.33 * 10 H r exp(-8.94r i/z ) , . IV-2,11 

where r la in pm and n(r) is in cm^yo" 1 . It ie. assumed that the phase function 

for A - 0.7 ym is valid at A - 0.8 ym. and that the albedo for single scattering iaj 

i 

S3 - .99999. The haze is assumed to have optical thickness values T - 0.5 and ; 
o 

T - 2.0. These optical thicknesses correspond to haze layers, which, if they are • 
M..0 km thick, are characterized by visibility ranges of ^3 km and n, 2 km respectively. 

The albedo of trees at A - 0.8 yo is A t ■ .53 and the albedo of soil is 

* j 

A « .05, see Krinov (1953). Four cases are illustrated in Fig. IV-2,2: i 

8 i 

_ I 

1. A small isolated bare soil patch exists in a forest so that A • ^-.53 I 

i 

A*- .05. Table 1 shows the radiances reflected from trees and j 

s i 

soil for several solar and viewing angles. | 

T 

2. A small group of trees exists in a region of bare earth so that A*A s *».05, s 

j 

A «• .53. Table 2 shows the reflected radiances as described above. j 

t 

3. A broken forest covers half a region of bare soil A m .29, “ .53, 

i 

A “ .05. Table 3 shows the reflected radiances. ; 

« 1 

4. A semi-infinite bare earth abutts against a semi-infinite deciduous j 

forest A “ .53, A - .05. The radiances emerging frcs tree covered j 

t a j 

areas are presented in tables la, and lc and those emerging fron bare 
soil areas are presented in tables 2a and 2c. 




> i 

























































































135 




:7:$V. 









Csaa 5 



Fig.IV-2,2 


Schematic views of various casia showing trees and soil 
which are considered in this study. 


Predicted radiances can be found by multiplying the tabulated values by 

F (the incident 6olar irradiance in the bandpass of the detector). For example, 
o 

if a spectrophotometer with a .01 Pm bandpass is utilized at A • 0.8 Mm, 

-2 ‘ -1 

Kondratyev (1969) gives the incident solar irradiance F ff “ 11.35 w a sterad . 


Examples of the downward directed diffuse irradiance functiou . 

1 

2 f D(A, t, p', p )p'dp', which appears in IV-2,6 are presented in table 4. 
> 0 ° 

IV-2»4 Discussion : 

Inspection of tables 1 through 3 reveal a number of characteristics s 
Comparison of radiances computed for P q ** .38 and P £ .38 in tables 1c 
and Id, 2c and 2d and 3c and 3d show that if the sun is near the horizon, 
ground features cannot be distinguished through a thick haze n«-"*r the horl~ 


I 


4 



I 

4. 

3 








zon even though such features may be evident near nadir. 

All of the tables show that for y ■ .38, the radiances reflected from 

o 

the haze at the solar azimuth angle are greater than those emerging at 180* 
from the solar azimuth angle. This is consistent with the observation that 
hazes produce a glare when viewed in the direction of rising or setting sun; 

it is ' a consequence of the forward peaked phase function of the haze. 

■ / 

/ * / 

All of the tables show that if y " .83, the radiances which are scat- 
1 o 

tered at <£>*,180° can be slightly larger than those reflected at the solar 

- ! ’ " - 

azimuth angle. The phase function used in these computations has a slight back- 

: Hr .! . 

scattering peak, see Deirmendjian (1969), which causes this phenomenon. 

!. ■ ' ' • • 

! : It is noteworthy that I - I , the difference between radiances 

I , K f 8 

reflected from trees and soil, for a given y aad y Q is found to be indepen- 
dent of azimuth angle. This is a consequence of the Lambertian character of 
these surfaces and it follows from IV-2,5 and IV-2,6. 


! - ‘*5 


Representative values oi the downward directed diffuse lrradlance function i 


which appears in IV-2,6 are presented in table 4. Values of this function i 

which are appropriate for various values of T, y Q and A can be obtained by inter- 

" 5 

polatlon of the numbers presented in that table. It should be noted that the 

earth's curvature renders results obtained for y < 0.3 somewhat unreliable. > 
, o 


: h 

: \ 


Visual contrast is defined by Middleton (1952) as 


A - 

i* 1 


C M “ " IV-2,12 

where 1 is the radiance reflected by an object and 1^ is the radiance re- ; 

fleeted by the background. The distinction between object and background is 

not clear if an inhomogeneous surface is viewed from above. Duntley et al. 

(1964) introduced a modulation contrast function which avoids the need to 
distinguish which surface feature comprises an object and which is the background: 


C R * **R “ *R^*R + I R ) 
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We will use C R In the following discussion. 

The contrast expected when the sun is high in the sky (y - .83) and the 

o 

ground is viewed near nadir (y ■ .96) has been computed to illustrate the effect 

of haze on a satellite image obtained near noon. This contrast, C , is plotted 

R 

against optical depth for the four cases considered in Fig. IV-2,3. Note that in 
the limit of ar optically thin haze, C_ * (A -A )/(A +A ). The small objects 
(cases 1, 2 and 3) are obscured much more rapidly by increasing the thickness 
of the haze than are the large objects considered in case 4. This is due to the 
veiling effect of the haze caused by light reflected from the surface, 
especially in case 1 where A » .53. 
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Fig. IV-2,3. Modulation contrast function, C , as s function of the optical 

K 

thickness of hazy atmosphere covering surfaces shown in FiglV*-2,2. 
The number on each curve identifies tbs csss considered. 
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Fig. IV-2,3 shows that the semi infinite representation of a horizontally 
inhomogeneous surface (case 4) overestimates the modulation contrast if the 
surface is actually- heterogeneous on a horizontal scale which is smaller than 
1 km. Multiple reflections between the hazy atmosphere and the surface degrade 
the modulation most severely if the surface is highly reflecting (see case 1 
with large A). 

It should be noted that the visual modulation contrast is proportional 
to the logarithmic response of the human eye; however , much data obtained 
from photoelectric remote sensing instrumentation is proportional to reflected 
radiances, I . Such data is obtained in digital .fora so that the algorithm 

fv 

which we have developed can be applied to digitized images even though it is 
aore sophisticated than logarithmic compression. 

IV-2.5 Conclusion: j 

* : 
A method to derive the total optical thickness of an atmospheric hare j 

i 

has been developed from radiative transfer theory. ‘ j 

If two recognizable objects whose albedos are known to be A^ and ora j 

found at zenith angle |i, if the mean albedo of the extended surface. A, can be 
estimated and if the solar zenith angle can be measured, then the contrast 
transmission function, ty* can be determined. If no prior knowledge regarding 
the albedos of surface features is available, but if some features are observed 
to be partially shaded from the direct solar radiation by a large object, then l|i can 
also be determined. The optical thickness which corresponds to the derived <j» can ba 
found from tab)e 4 and IV-2,6. Measuring the reflected radiances at various 
wavelengths will therefore yield the optical thickness as a function of the 

j 

corresponding wavelengths. j 


/• ..j 


/ 



/ - 



We wish to thank Dr. J. K. Hansen of the Institute for Space Studies G.b.F.C. 


for providing us with a. copy of his doubling program. 


References for Appendix IV-2. 


Breitling, P.J., PilipovskyJ, S.,(1970): Computer Simulation of Optical Con- 
. j trast Reduction Caused by Atmospheric Aerosol. AIAA paper #70-194, 

, AIAA 8th Aerospace Sciences Meeting, New York. American Inst, of 
Aeronautics and Astronautics, 1290 Ave. of Americas, N.Y. 

i " I 

Deirmendjian,(l969) : Electromagnetic Scattering on Spherical Polydlsperslona, 
American Elsevier Publishing Co., Inc., New York, pp. 290. 

I I 

Duntley, S.Q., Boileau, A.R. , Prelsendorfer, R.W., (1957): Image Transmission 

by the Troposphere I,J. Opt. Soc. Aa. 47 , 499-506. 

Duntley, S.Q., Gordon, J.I., Taylor, J.H., White, C.T., Boileau, A.R., 

Tyler, J.H., Austin, R.W., Harris, J.L. (1964): Visibility App. Opt. J3, 

549-598. 

\ 

Hansen, J.E.,(1969): Radiative Transfer by Doubling Very Thin Layers, 

Aetrophys. J. 155 , 566. 

Irvine, W.M., 0.968): Multiple Scattering by Large Particles., Aatrophys. J. 

152. 824. ... . . , 

Krinov, E.L. ,(1953): Spectral Reflectance Properties of Natural Formations, 
translation TT-439, Natl. Res. Council Canada, Tech., pp. 268. 

Kondratyev, K. Ya., (1969): Radiation in the Atmosphere . Academic Press, Inc., 

New York, App. 4, pp. 912. 

Kondratyev, K. Ya., Buzinkov, A. A. , Vasilyev, O.B. and Smokcy, 0,1. ,(1973): 
Influence of the Atmosphere on Spectral Radiance and Contrasts of 
Natural Formations Measured from Space. Proceedings of U.C.L.A. 
International Conference on Radiation and Remote Probing of the 
Atmosphere . J.G. Kurlyan, Editor. U.C.L.A., Los Angeles, CA, pp. 309-336. 

Middleton, W.E.K. (1952): Vision Through the Atmosphere , University of 

Toronto Press, Toronto, Ont. pp. 60-64. 

Van de Hulst, 0971): Some Problems of Amlsotropic Scattering in Planetary 
Atmospheres, presented at Symposium #40 on Planetcry Atmospheres, 

Marfa, Texas, Oct. 26-31, 1969. C. Sagan, T. Owen, H. Smith, Editors. 

D. Re id el Publishing Co., Dordrecht, Holland, pp. 177-185. 


140 


j~ \ •» - 


Appendix IV-3: The Direct Use of a Nitrogen Laser for a Broad Sand High 

Resolution Lidar (System IV) 

IV- 3a) Experimental Considerations of System IV 

In preliminary experiments we showed that individual lines from the 
laser output at 3370 X can be as narrow as about 4mX; however the energy 
is distributed among a large number of 1 ines covering an approximately 1 X 
bandwidth. Consequently, using conventional practices, an impractically 
small fraction of the Nj laser light is available for the high-resolution 
scheme. 

It would be possible to make use of the full output energy if the output 
could be directed to obey the prooer wavelength vs. angle law. An explanation 
and demonstration of this scheme is contained in the papers of Trauger and 
Roesler (1972), Roesler and Stoner (1973), and Trauger et al.(1973). 

Basically the scheme is to match the change in wavelength with viewing angle 
to the angular dispersion of the spectrometer. As shown in the references, 
excellent results are obtained when the wavelength vs. angle is produced by 
source motion. To apply this idea to the atmosphere, it is necessary to 
illuminate the atmosphere with light having the proper X vs. 0 behavior. 

What needs to be demonstrated is that one can force the laser output 
to deliver its energy In the proper A vo 0 relationship. Since a special 
laser would have to be built for this purpose, we have investigated the 
possibilities using an pumped dye ceil. The A vs 0 function was obtained 
using a wedge etalon as shown in Fig. IV-3,1. 
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Fig. IV-3,1. Schematic view of a system IV transmitter. 

Light passing through the top of the etalon has wavelength X 2 * reflects 
from the mirror back through the same part of the etalon and is directed 


by the lens through the active region of the dye cell at angle Similarly 

light of wavelength A^ passes through the lower portion of the dye cell and 
finally emerges at angle 6^ ; intermediate wavelengths are found at inter- 
mediate angles. Figure IV-3,2 shows a photograph of the output light dis- 
persed horizontally by a coarse grating. Five grating diffraction orders 



Fig. IV-3,2. Output of laser system IV prototype. Output angle is 
displayed vertically, while grating dispersion is horizontal. Dis- 
persion in nighest (right-most) order is 'v,5X/mm. Lower exposure is 
for a smaller wedge angle. 

of the active range of the dye are seen, each crossed by several diagonal 
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been changed between exposures, and consequently ;he wavelength difference 
between top and bottom of the etalon is different.' While this was a low 
resolution experiment for demonstration purposes only, the application of 
the principle to high resolution is straightforward. Fig. IV-3,3 shows 
the set-up used to obtain the results. 


rn • p* UBfc lOdt* _.T3 t=S.:^^!T; 




Fig. IV-3,3. View of N^-pumped dye laser system configured as 
system IV type transmitter. Note output projected on wall. 


To make this scheme work with the laser one must use a slightly 
different approach, the elements of which have been demonstrated. As shown 
above, to achieve the proper X vs.0 behavior, the active medium must be 
directed. Since the dye cell is short and can amplify at any X within its 
range, an etalon is sufficient for this. The laser is highly super- 
radiant and photons generated spontaneously at one end are amplified as 
they travel down the discharge tube. In order to prevent this, properly 
selected "leader photons" must be present in the laser tube before it 
becomes active. This can be accomplished by the following proposed 
oscillator-amplifier scheme: Two lasers are triggered by the same spark 
gap trigger, the first laser being adjusted to become active 3-4 ns before 
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the second. The output of the first is filtered by an external passive 
etalon to obey the proper X vs. 0 function and fills the second with leader 
photons by the time it becomes active. Thus a high output is obtained 
with about half the overall efficiency of a single, undirected laser output. 
The oscillator-passive etalon-amplifier has been satisfactorily tested 
using nitrogen pumped dye cells as shown in Fig. IV-3,4. 


X vs. 0 
output 



output 


beam-splitter 


Fig. IV-3,4. Simplified schematic view of oscillator-amplifier dye 
laser for system IV use. 
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The experiments with the N 2 laser-pumped dye cells have been conducted 
to demonstrate the principles and gain experience with the various possi- 
bilities before designing a special N 2 laser system. The experiment 
illustrated in Fig. IV-3,4 has been conducted with a large gap etalon 
providing high resolution. We have also conducted experiments to determine 
the optimum timing of pulses in the dye laser - amplifier combinations. 

Since the laser pulse is short. It is important that the pumping radiation 
reach the amplifier dye cell properly timed so that maximum efficiency is 
obtained in amplifying the narrowed etalon output. The delays caused in 
narrowing the oscillator output by the grating-etalon combination are 
significant,, and the scheme shown in Fig. IV-3,4 is not optimum; greater 


1 J' 








amplifier dye cell. 

The status of experiments aimed at optimizing the timing is as 
follows: we have noted the delays and the loss of efficiency when the 
separation between the grating and dye cell is increased, and have made 

1 ■ 

' some attempts to increase the efficiency by adjusting the timing. ->e 

efficiency experiments have been somewhat inconclusive. Partly there has 

I I ‘ 

I been insufficient time to do the careful experiments necessary and partly 

ji • 

I there is. a need to make several new mountings, obtain new optics with 

I i i ! 

proper specifications, and make modifications to the dye laser apparatus. 




These experiments will continue, and progress may be expected in the 
coming months. 

Experiments are also being made to determine the optimum dye cell 
size and the optimum splitting of the Nj laser light. Considerable progress 
has been made, but more careful work under more carefully controlled 
conditions is required. 

IV- 3b) Theoretical Limitations of System IV . 

The broad band high spectral resolution lidar system is designed so 
that the radiance which is propagated at an angle with respect to the 
laser axis is a unique function of the wave length. Similarly the radiance 
which is backscattered to the lidar receiver is expected to be a unique 
function of t|>', the angle between the propagated direction and the receiver 
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where c is the wave number. 

The resonance condition in an etalon with spacing £ and refractive 
index li requires 

2y£o cos0 ■ k 

where k is an integer, and 6 is the angle between the optic axis and the 
direction of propagation. For two adjacent angular zones n and n+1. 
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a cosfl “ o . ,cos0 ,, 
n n n+1 n+1 

which can be rewritten, for 8« 1, 

‘I ■ a n 9 n °n+l 0 n+l 
°n ‘ °n+l ' ~2 2 


26a 2, „ fl2 n 2 

~ R - 6 n+l 0 n 


IV-3,3 


0 » 0 , 0. » 7.7 x 10 rad 
o 1 
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0 ** Jn 0 

n 1 


If the etalon system is designed to accept 20 zones, i.e. n * 20, then the 


maximum receiver field of view angle ■'.3 

9 20 " 3 ** 5 x 10 3 radians 


lV-3,5 


The ratio of apertures of the telescope and the etalon reduce this angle 

so that for a 5 cm diaiu. etalon system and aim diam. telescope which is 
proposed for the Shuttle 

'P' « 1.72 x 10~ 4 radian IV-3,6a 

while an aircraft system with a 20 cm diameter telescope has 

i}*' - 6.90 x 10 _ * radians. IV-3, 6b 

Conditions cited in Eq. IV-3,1 can be achieved if single scattering is 
the only scattering mechanism; however, multiple scattering degrades the 
angular and hence the spectral resolution of the proposed high spectral • 
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resolution lidar system. 

The theory developed by Weinman (1975) was used to estimate the magnitude 

of multiple scattering. The hazes which ve Investigated were assumed to have 

extinction coefficients $ ■ 2 km * and 0.2 km ^ which are independent of 

altitude. These models correspond to visibilities of 2 and 20 km respectively. 

These conditions are found at the surface for the BF - 10 and BF * 1 models 

respectively; see Tables 3-1 and 3-2. The Ildar is assumed to be located 

distances h ^ from the haze; h ■ 185 km for the Shuttle case and h * 10 km 
o o o 

for the aircraft case. The phase function is approximated by the parameters 
applicable to Haze-C cited in Table 1 of the Weinman paper. and are 
the total and the singly scattered lidar return measurements. 

Table IV-3,1. Ratio of Multiple to Single Scattering, , from 

a System IV Lidar. 



6=2 km -1 

6 = 0.2 km -1 

Shuttle (b = 185 km) 
o 

.213 

.030 

Aircraft (h = 10 km) 

..065 

.008 

o 




(1) The distance separating the- lidar from the haze, b , is designated H’ 
by Weinman. 
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The entries in Table IV-3,1 represent the proportion of multiply scattered 
lidar photons which both enter the atmosphere within the angle if/ * (Eqns. IV-3,6) 
and return within the same angle. Since any multiply scattered photon has 
roughly a .95 probability of being received in a different resolution zone, 
this has the effect of degrading the high spectral resolution and nullifying 
the potential advantage of this technique. 
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In summary, the scheme presented here for producing a wavelength encoded 
angular distribution of laser radiation, together with a properly tuned 
PEPSIOS receiver, make it possible to use a larger proportion of the 
available direct Nj radiation in a proposed HSRL system. In the example 
presented here, 2X of the available energy available across a 1 X band 
* width, has been utilized at an effective resolution of 10 ^ X. However, 
the solid angle of the field of view has been correspondingly increased 20- 

i ' 

fold. Thus solar background radiation has also increased 20-fold. This 

j ' \ ‘ ‘ : 

f scheme | is susceptible to degradation of the spectral resolution due to 

I ' r 

! i 

multiple scattering within the field of view, as discussed above. Finally, 
the tolerance in angular alignment of the laser and receiver field of view 
has been decreased by a factor of 1/2/20. 
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5-0 Conclusion 

The most widely dispersed aerosol in the atmosphere originates over 
deserts. Models of the desert aerosol were defined on the basis of the 
varied data presented in the literature. Radiative heating and cooling rateu 
produced by such model aerosols were computed using computational schemes 
which are sufficiently fast yet accurate enough to be Incorporated into 
General Circulation Models. (The results obtained from General Circulation 
Models which Include perturbations due to desert aerosols are not presented 
in this report. The first results are emerging as this conclusion is being 
written.) General Circulation Models are Inherently noisy due' to the finite 
temporal and spatial increments which they employ. Such noise corresponds 
to an uncertainty of 0.5°C/day in the heating rate. An analysis was con- 
ducted to show how this uncertainty in heating rate ultimately manifests 
Itself on the aerosol optical extinction profile. This analysis provided 
the foundation for the design of two specific lidar systems: one which could 

be flown on the Shuttle missions and the other on an aircraft-. 

High Spectral Resolution Lidar (err HSRL) systems may be used to sep- 
arate the light scattered by aerosols from that scattered by molecules in 
order to determine vertical profiles of the aerosol extinction coefficient. 
All HSRL systems considered were eyesafe for the unaided eye looking up from 
the Earth's surface. 

It was found that the HSRL dye laser experiment is feasible for air- 
craft deployment using present technology. Solar noise is not a problem due 
to the relatively high returned powers, and consequently HSRL operation in a 
Fraunhofer line is not necessary. 
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The simplest aircraft borne system would deliver 3 mw of optical 


power @ 60 Hz from a platform located at h^ " 10 km above the Earth's surface 


Such a system will provide values for the aerosol extinction coefficient in 


vertical Intervals of depth u ■ 0.75 km with sufficient accuracy to determine 

-1 


atmospheric heating due to aerosols within ± 0.5°K day 


i ■ 


The simplest Shuttle borne system would deliver 2 dm of optical power 


@ 60 Hz from a platform located at h Q “ 185 km above the Earth's surface. 


! This system also determines the aerosol extinction coefficient in vertical 


fi 


intervals of depth u “0.75 km with sufficient accuracy tc determine atmo- 

-1 


spheric heating due to aerosols within ± 0.5°K day ~. This system yields 
solar background levels which are higher than the returned signal strength; 


however operation in a Fraunhofer line using a receiver field of view which 

is a practical optical (component) limit gives a feasible value for N /N . 

s m 

Thus, the Shuttle borne HSRL is theoretically capable of providing the 
desired information on , but it is at the limit of current technology. 

Solar noise could be reduced by employing a flash-lamp pumped 
transmitter which would deliver 3 mw @ 1 Hz repitition rate with the desired 
spectral resolution. This would eliminate the need for a very s.iall receiver 
field of view and operation in a Fraunhofer line. However it would be diffi- 
cult to attain both the high spectral resolution and high spectral stability 
using a flash-lamp pumped transmitter. 

It appears that a 15 mv @ 60 Hz dye laser is forthcoming. Our 
analysis shows that HSRL operation from the Shuttle using such a dye laser 
system, operating in a Fraunhofer line, is capable of yielding vertical 
profiles of the aerosol extinction coefficients both by day and night 





